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Preface 

It has been noted, in the course of study of geodesy, that 
there is a sparsity of material devoted to the mechanics of determining 
the size and shape of the reference ellipsoid from actual observations. 
Bare outlines of some methods are given in the few current textbooks 
available in English, and although some of the older texts, notably 
Clarke (5) and Crandall (6), do give examples, the methods used are 
rather ancient and outmoded and in notation systems quite unfamiliar 
to the present-day student. Hayford’s original study of 1909 is the 
best available (10 ). Since one of the basic aims of geodesy is to 
determine the figxire of the earth, it seemed like a good idea to 
pursue the matter, if even yet somewhat superficially, to the point 
where there would be few questions left concerning the manner of going 
about the task. 

The writer wishes to acknowledge the help and advice provided 
him in preparation of this thesis, particularly by his adviser. 

Dr. W.A. Heiskanen, Director of the Institute of Geodesy, Photo- 
grammetry, and Cartography, the Ohio State University. Acknowledgement 
is also due Dr. U.A. Uotila and Mr. Ivan Mueller, both of the Ohio 
State University, for certain helpful sviggestionsj the U.S. Army Map 
Service and Mr. W.M. Kaula, for data on points in the Finland - 
South Africa arc| the U.S. Navy, without the instigation of which, 
the course of study referred to in the opening sentence would not have 
been possible. 
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Finally, the writer wishes to express his appreciation to 
Frauke and Peter Wilson, the former for her fine job of typing the 
bulk of this thesis on short notice, and the latter for certain 
editorial assistance to the typist. 
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Chapter 1. Introduction 



1.1 Preleminary Remarks. 

That the most convenient practical reference surface for purpose 
of geodetic computation is the ellipsoid of revolution has long been 
recognized. This is emphasized by the fact that the major datums in 
use in the world today are all based on ellipsoids of revolution. The 
fact that this shape is a relatively simple regular mathematical sur- 
face, and that it closely approximates the mean geoidal form, make its 
choice quite understandable. The Russians, who have been proponents of 
a triaxial ellipsoid, do not even consider its use for practical 
geodetic computations. Izotov (l8) states 

"The trieixial ellipsoid could not be accepted for triangulation 
processing, as this would have made the working out of geodetic 
measurements too complicated". 

To dismiss the subject of triaxiality for purposes of this paper, the 
following is cited from "The Earth and its Gravity Field", by Heiskanen 
and Vening Meineszs- 

"It appears that we had better forget the triaxial ellipsoid, 
particularly because since the earth is in close isostatic 
equilibritim, it can hardly have the shape of a triaxial 
ellipsoid." (13) p. 79"8o. 

The use of geodetic data to determine relatively precisely what 
size and how much flattened the reference ellipsoid should be, is the 
subject of this discussion. In this connection, one of the principal 
purposes of geodesy has been variously stated; - 
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"Its theoretical function is to determine the size and shape 
of the earth ...” (I3)j "Geodesy is the science which treats 
of investigations of the form and dimensions of the earth's 
surface." (17) I "The accurate determination of the figure of 
the earth.”' (6); "In combination with observations for latitude, 
longitude, and gravity to assist in determining the size and 
shape of the earth, ..." (2); "It may be said that the ultimate 
aim of scientific geodesy is to determine the size and shape of 
the geoid, ..." (2). 

Of course, data collected in geodetic surveys along with the calcula- 
tions connected therewith, and computations and results of computations 
of the figure of the earth al*e mutually interdependent. 

Although the problem of determining the best size and shape 
for the reference ellipsoid has decreased in importance relative to 
other modern forms of geodetic endeavor, there are still significant 
efforts being made in this field, notably by the U. S. Army Map Service. 
This is not to say that the exact dimensions of the earth are not im- 
portant, but the position has now been reached where further precision 
given to a reference ellipsoid must be subordinated to the problem of 
determining the deviations of the actual form of the geoid with respect 
to the reference surface selected. 

This change of emphasis is described somewhat succinctly by 
Tengstrom; - 



"The principal problem of geodesy is to determine the geocentric 
coordinates of all points on the earth's surface, a task which 
was earlier referred to as a determination of the 'size and 
shape' of the earth." (30)j 



or by L. G. Simmons s- 



"The ultimate goal of the geodesist is, I suppose, to determine 
the parameters of an ellipsoid of revolution which best fits the 
figure of the earth as a whole. But this is not all. He also 
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is concerned with the details of lack of fit of this ellipsoid 
to the actual Earth's shape, the geoid." (28). 

Most important would be the selection of one reference surface 
to be used in problems of world -wide extent, such as for requirements 
of missile warfare, the success of which might very well be determined 
by accurate knowledge of locations of launch sites and targets. Natu- 
rally, the best reference surface obtainable is to be preferred. The 
establishment of a "World Geodetic System" and conversion of existing 
systems to it, has for several years been one of several significant 
projects promoted and encouraged by Dr. Heiskanen (l6) (13) p. 299-310. 
Both the U. S. Air Force and the U. S. Army Map Service are developing 
world geodetic systems for military purposes. These are eventually to 
be combined into a single unified Department of Defense World Datum, 
the contributions to the establishment of which will be made by all 
U. S. military services. The reference ellipsoid to be used for this 
will be based on all available information from all parts of the world. 
( 3 ) P. 77 . 

Having exposed the near-obsolescent character of the problem 
of determining the dimensions of the reference ellipsoid, we may now 
proceed to describe the several variations of its solution. 

1.2 General Discussion of Methods. 

The standard approach to the problem under consideration has 
been, from earliest times right to the present day, that contained in 
the arc method. Variations of this approach have become refined and 
quite sophisticated, taking into account the earth's isostatic 
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equilibrivmi, and in some cases being combined with the gravimetric 
method. The earliest recorded application of the arc method is told 
in the now trite story of Eratosthenes* activities in this field. 

The arc method assumes that short meridional arcs on the 
earth's surface are arcs of, circles; when two meridian arcs are 
measured, one at relatively high latitude, and one near the equator, 
both the equatorial radius and the flattening can be computed. The 
more such arcs used, combined in a least squares solution, the better 
accuracy is obtained. 

A more convenient adaptation of the arc method has been used 
since early in the present century. This uses a cmbination of 
geodetic and astronomic data in the form of astro-geodetic deflections 
of the vertical. Hayford used this method in developing what is now 
known as the International Ellipsoid (lO). He called the procedure he 
used the "area method", because the data was taken from many points in 
a large network of triangulation. He made several solutions, the 
deflections of the vertical in each being reduced to the geoid (co- 
geoid) according to the theory of isostasy using different depths of 
compensation. He chose the solution in which the svun of the squares 
of the residual deflections after adjustment was the least. Thus a 
"by-product" of his work was a solution of the mean depth of compen- 
sation for the earth's crust. 

The application of the gravimetric method to the problem under 
consideration is really a combination of the arc method, modified by 
the introduction of gravimetrically determined deflections of the 
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vertical, to obtain the earth's equatorial radius, and the deter- 
mination of the earth's flattening from gravity anomalies by use of 
Clairaut's formula. 

The latest, and of course the most modern method used to 
determine the earth's shape, is the artificial earth satellite method. 
5y measuring the perturbations in the satellite's orbit around the 
earth, caused by deviations from the spherical in the earth's shape, 
the flattening is obtained. , This method is so new in actual 
application, that an evaluation of results thus far obtained can only 
be the subject of conjecture. 

One other method, similar to the deflection of the vertical 
application of the arc method, involves the use of 'the geoid un- 
dulations, obtained both graviraetrically and astro -geode tically. The 
solution is obtained in reducing the sum of the squares of the 
differences between the two at many points to a minimum. 



Chapter 2. Dimensions fiom Two Arcs 



2.1 Meridional Arcs. 

In this analysis, as in all subsequent ones in this paper, the 
reference surface is assumed to be an ellipsoid of revolution. The 
mathematical properties of an ellipse fom the basis for expressions 
referring to the meridian of such a surface. 

The measurement of two arcs in meridian, one near the pole and 
the other near the equator, provides the data for this, the simplest 
means whereby the earth dimensions may be obtained. The arcs, which 
may be on different meridians, are measured, and the difference in 
latitudes between extremities for each is computed, the latitude of 
one end point for each arc having been determined astronomically. 




In fig. 1, <p^ and cp^ are the mid-latitudes, and s^^ smd s^ the 
measured lengths of the two arcs. Assuming s^ and s^ to be arcs of 
circles, we can write s- 
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( 2 - 1 ) 



R, = T— CSC 1" 
1 



«2= 



and Ai:p^ being expressed in seconds. 

The expressions for the radii of curvature of the elliptical 
meridian at the respective latitudes are : 



( 2 - 2 ) 



R-. = 



- ^(1-0 



1 " 2 . 2 x3/2 

(1-e sin cp, ) 



^2 = 



a(l-e^) 



(l-e^sin^cp^)^'^^ 



(2-3) 



Equating the two expressions (2-1) and (2-2), we obtain: 

s. 






esc 1" = 



a(l-e^) 



(1-e^sin^cp^)^'^'^ 






CSC 1" = 



a(l-e^) 



/-, 2.2 v3/2 

(1-e sin cp^) 



2 

Dividing R^ by R^ in the form (2-3) and solving for e‘ 



(2-M 



Sj_^2 



= <1 



2 

e = 



2 2 

1-e sin cp^ 

2 2 
1-e sin cpj^ 



2 2 2 
sin cp 2 "'l 



Then, the value of a can be found from (2-3) : 



(2-5) 



a = 



,, 2 . 2 
s^(l-e sin cp^) 

Acpj^(l-e^)sin 1" 



S2(l-e^sin^cp2)^''^ 

Acp2(l-e^)sin 1" 
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Now all that remains is to determine the flattening, which is 
simply: 




2.2 Oblique Arcs. 

In the early efforts of measuring the earth's figure relatively 
short arcs were measured, almost always in meridian. In later deter- 
minations, existing triangulation data was put into use, and consequent- 
ly successive points used were not always on the same meridian. Such 
geodetically computed arcs, chosen so as to be inclined to the meridian 
at as small an angle as possible, could be reduced to the meridian by 
appropriate formulae or from tables based on the used ellipsoid. Thus 
armed with the same data as in section 2.1 a similar problem can be 
solved. 

When the length of the oblique arc is measured, along with its 
azimuth and mid-latitude, its meridian distance is given by: 

12 2 

(2-7) AcpR sin 1" = - s cos A - S sin A tan cp 

+ sin^A cos A (l+3tan^cp) 

6ir 

where W is the radius of curvature of the prime vertical and A is the 
azimuth, reckoned positive clockwise from south. 

Using (2-7) a chain of triang!.es can be projected to the meridian 
for the whole length of the chain. Given geodetic latitudes, however, 
use can be made of specially prepared table, such as (3^)^ (35)^ or 
( 36 ), which have meridian distance tabulated against latitude for the 



ellipsoid involved. 
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The above method of detei’mining the earth's dimensions (Sections 
2.1 and 2 , 2 ) ^ though it represents the basis for other methods, is 
mainly of historical and mathematical interest, since no longer are 
serious scientific computations carried out in such simple terms. 
References for the presentation given in these two sections are (6), 
(17), and (23). 



Chapter 3- Dimensions from Several Arcs. 



3.1 General. 

When geodesists of the l8th Century first began to compare the 
results of their arc measuring efforts, they found such large discrep- 
ancies that they began to speculate about the accuracy of their meas- 
urements and about the correctness of the assumption that the earth's 
tme shape is an ellipsoid of revolution. Of course the accuracy of 
surveys of those early days was not what it is today. Nor could one 
hope to obtain the same results from arcs, taken two by two in dif- 
ferent parts of the world having varying kinds of topography, meas- 
ured by the crude methods of the day. Not only was the measurement of 
geodetic triangles so rude that spherical excess remained xindetected, 
but the process of adjustment by the method of least squares was un- 
known (23)- 

Obviously, not wanting to discard such a convenient reference 
surface as an e'llipsoid, which after all departs from the geoid sur- 
face by probably less than 100 m. at the maximum, and from the 
equilibrium figure by even less (of the order of 2 m. at cp = 0,45, 
and 90, (13) P* 59), the best thing to do would be to combine the 
results of several arc -measurements in such a way as to determine the 
mean values for the earth as a whole. 

3.2 Laplace's Method. 

Toward the end of the l8th Century, Laplace published the re- 
sults of his attempt at combining the data from nine different sources, 
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extending from Lapleind to Cape of Good Hope in latitude. His method 
involved the data corresponding to an arc of one degree of latitude for 
each of nine different meridians. 

Again considering these one degree arcs as arcs of circles and 
using R, the radius of curvature of the midpoint of the arc, and cp, 
the mid-latitude, we have:- 



( 3 - 1 ) 



d = 



2 T/ R ^ TT a (1-e^) 

l8o (1-e^ sin^ 



Developing this in series, we obtain; - 

/o ^ Tra(l-e^) . 3 2 . 2 15 4 , 4 s 

(3-2) d = (1 + — e sin cp + e sin cp + . . . ) 

which can be written in the following form;- 

( 3 - 3 ) 

M 



2 4 

d = M + N sin q) + P sin 9 + ... 



in which 



- 'na (1-e ) 

iBo 



3 2 

N = I e M 



15 4 

P “ ^ e M, ... etc. 



Powers of e higher than the square are neglected. 

In order to combine the data for each degree -measurement in 
the manner yielding the least discrepancies, the observation equations 
for each arc were of the form;- 

( 3 - 4 ) d^ - M - N sin^ 9 = v^ 

2 

dg - M - N sin 9 = v^ etc. 

in which v^, v^, etc., are the residual errors. Since the method of 
least squares was unknown at the time, Laplace combined the error 
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equations in such a way that the absolute values of the residual errors 

would be a minimum, and the algebraic sm equal to zero. Then solving 

2 

simultaneous equations for M and N, the values of e and a are found. 

Of course, this problem is easily adapted to least squares solution. 

The results of Laplace's computation aret- 

M = 11 06i. 2 meters 
N = 1 19^.9 meters 

e^ = 0.007202 

a = 6 383 600 meters 

f = 1/278 

The error, or difference between a measured degree and one 
computed from the result of the crude adjustment of the nine arcs, is 
of the order of 268 meters. From this, Laplace concluded that the 
earth deviates materially from an elliptical figure ( 23 ). 

3«3 5 y Correction of Latitude. 

A refinement of Laplace's method, in use by the end of the 19th 
Century, is here described ( 6 ), ( 17 )* This also involves the formation 
of observation equations between observed latitudes and either directly 
measured meridian arcs or oblique arcs projected to the meridian. 

For short arcs the relation s = AT R sin 1 " can be used, but 
for longer arcs a correction is required. Crandall ( 6 ) states that the 
simple relation given is good for a Acp of "several degrees" because of 
the large probable error of latitude determination, some oVo 4 (more 
likely this will be on the order of 0"1 at best), in addition to the 





il 



13 



error caused by local deflections of the vertical. The correction, for 
use with longer arcs is:- 

(3-5) d s = ^ a e^ (/!yp" sin 1")^ cos 2 cp 

Given astronomic latitudes for several stations in different 
arcs, as well as geodetically measured segment lengths £, we choose 
separate initial points for each arc. From the first two latitudes in 
any one arc we have;- 

(3-6) cpg - cp^ = sin - R sin l" 



where 

(3-7) 



1 ^ (1 - e^ sin^ cp)^^^ 

^ a (1 - e^) 

In this method we make use of the used reference ellipsoid, for 



which we know the values a^ and e^ , and. compute corrections 6a and 

6e , to that ellipsoid resulting from the data given. Thus a = a^+ 6a, 

^2 2^2 
and e = e + 6e . 
o 

Expanding equation (3-7) into a series and substituting it in 
equation (3-5) t- 



(3-8) cpg - cPi - sin 1" 



1 ^ 3 . 2, V 6e 

R -2 + (1 - 2 'P) — 



a 



ds 



R sin 1" 
o 



in which R is the value of R for a and e . 

o o o 

Considering the measured meridional arcs accurate in comparison 



with the observed latitudes, we give corrections v^, v^, to latitudes 

r 

cPl, cpg, etc. Since in the derivation of (3-8), all terms containing ^ 

2 2 

were neglected in the coefficients of 6 a and 6 e , i.e., e was con- 
sidered equal to zero, we can say R = a and s/(a sin 1") = cp^ - cp^. 
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Thus equation (3-8) becomes:- 

(V'a) - = R sin 1" * a sin l" [’ 



d s 



R sin 1 
o 



TT ) and then 



(3-9) Vg = 



Tg - 9. 



a sin 
o 



^ 6 a + (l - I sin^ cp) (cpg - (f>^) 8 e‘ 



ds 



R sin 
o 



^ - (^2 ■ ^Pj)’ 



This becomes the basic observation equation to be written for 

2 

each segment of each arc. The ttnknowns are v^, 5 and 6 e . 

2 

In the least squares solution, in which S v = minimum, the 
magnitudes of the various coefficient elements make it convenient to 
use the following modifications. 

p 

Placing X = 6 a/ 1000 y = 1000 5 e 

m^ = - 1000 (cpg - 9^)/% 



92 - 9i 
1000 

s - ds 



(1 - I sin cp) 



^2 R sin 1" " ^*^2 " “^l^ 

o 

We can then write a simple error equation for each segment, 
substituting the above symbols in (3-9) s- 



Vg . + ^2 X + ng y + Ig. 

Since the initial latitude must be corrected, an additional 
observation equation is written for it, simply v^ = v^. In this 
method, there will be as many unknowns as there are arcs plus two, i.e., 
the V for the initial latitude of each arc, plus x and y. In the case 
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that a value for the initial v is assumed to he some arbitrary value, 
e.g. zero, there will be only two unknowns. 

For each arc, observation equations will be as follow;- 



( 3 - 10 ) 

-nn^x + jyr + = Vj 

Yl + mjX + njy +13=^3 
* * * * * 



v^ + m X + n’y +1 = v 

In n*^ n n 



The result of a computation of this type presented by Crandall 
(6), and using data (circa 1890) given in the "Handbuch der Vermessungs- 
kunde" ( 37 ) for six arcs (French, English, Hanover, Prussian, Russian, 
and Swedish ) , varying in length from 0?5 to 12 ? 4 , is as follows ; - 



a = 6 377 800 meters 

e^ 0.0069091 
f = 1/289 



The standard error of the mean amounts to 2 Vl for a latitude 
determination, which is quite a bit greater than the accuracy obtain- 
able for first order work, (about O'.'l). Thus, Crandall concludes ( 6 ), 
an ellipsoid of revolution will not fit the data without large local 
deviations of the vertical. 



3 »^ By Correction of Radius. 

A method given by Bomford (2) p. 364 , which is a variation of 
that described in Section 3 * 3 , works with radii of curvature at points 



in an arc rather than with latitudes. 



l6 



For given values of R in widely separated latitudes, the values 
2 

of a and e are obtained from the expansion of equation (2-2) 



(3-11) R = 



a (1 - e^) 
(l « e^ sin^ 



= a (l - e^) (l + g- e^ + ^ e^ sin^ cp + . . . ) 

2 2 2 
Using e = 2 f - f and the approximation 6 e « 2 5 f, and 

2 

dropping the terms of higher order than a , the following form is 
obtained;- 



( 3 - 12 ) 



6R=6a-2a &f(l-J sin^cp) 
O d 



in which 5 R = R - R results in changes 6 a and 6 f in a and f f 
the el-lipsoid provisionally adopted. From equation (2-3) we may use 
R = s CSC 1”/^" for arcs shorter than about 100 miles (2). The 
value of R^ would be obtained for the used ellipsoid. 

In a more practical application, more than two values of R 
would be available, and observation equations for least squares 
solution would be as follow:- 

S CSC 1** 

(3-13) v^ = X - 2a^y (l ” g sin cp^) ^cp^" 

, S CSC 1" 

v_ = X - 2a y (l - § sin cp_) — — n — + (R„) 



ACp. 



0^2 



***** 



V = X 

n 



^ p S CSC 1" 

2a y (l - ^ sin'^ cp ) — + (R ) 

o*' ' 2 ^n'^ Acp ' o'^n 






where x = 6 a and y = 6 f are the unknowns, s^ is the measured length 

2 

of the arc Acp^" = cp^ " seconds, etc. Again Z v is 

made a minimum. Weights may be assigned as deemed necessary. 



Chapter 4. Dimensions from Arcs, Using Deflections 

of the Vertical 

4.1 Basic Considerations - Astro-Geodetic Deflections. 

(Note; A more complete analysis of "basic considerations," in- 
volving absolute deflections of the vertical, is given in section 4.3.) 

As mentioned earlier, a very convenient adaptation of the arc 
method makes use of astro-geodetic deflections of the vertical, which 
are generally available in profusion wherever triangulation control nets 
have been made. This has been done in many parts of the world in 
connection with the various national surveys being carried out. 

This method can be applied to individual arcs of meridians in 
one extreme, and to the whole network of triangulation points in the 
area covered by a national control survey in the other. Because of the 
latter application of the method, Hayford chose to ceill it the "area 
method" when the results of the U.S.C. & G.S. computations of the figure 
of the earth based on data collected in the large area covered by the 
United States were published (lO). 

Here again, the reference ellipsoid used in computation of the 
geodetic data is used as a point of departure, or more accurately, an 
"area of departure", from which we must apply measurements and compu- 
tations to get as close to the best ellipsoid as we can. This "best" 
ellipsoid is the one most nearly approximating the earth spheroid or 
mean geoid surface, in size, shape, and orientation. 

Figure 2 shows in simple terms, the relationship which catn exist 
between a reference ellipsoid and the geoid. 
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or- 




Fig. 2 - Astro -Geodetic Deflections of the Vertical 

In general, there are six different items that may be corrected 
by the procedure involved, and the particular ones that are corrected 
depend upon what is desired and upon what and how much data is supplied. 
The existing geodetic data, i.e., latitudes, longitudes, deflections, 
etc., depend upon the initial point (datum) of the survey, and upon the 
assumptions made about that point. The six items are: the equatorial 

radius aj the polar flattening f; and the following quantities assiAmed 
for the initial point: its latitude cp^, its longitude the geoid 
height N^, and the azimuth A of its meridian. In many cases, all that 
are desired would be the corrections Sa and Sf. Usually the cor- 
rections to cp and \ at the initial point will be in the form of 
^o o 



corrections to the meridian and prime vertical components to the 
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deflections of the 
The following 
cpS X' 
cp, X 

(cp^-cp) 

(X^ -X)coscp 

f 

V 

dC 



D 



< 



vertical. 

notation will he used in the discussion to follow: 

- Astronomic latitude and longitude 

- Geodetic latitude and longitude 

- Meridian component, astro-geod. defl. 

- Prime vertical comp., astro-geod. defl. 

- Residual deflection in meridian 

- Residual deflection in prime vertical 

- Orientation corrections to deflections 

due to displacement of ellipsoid 

- Change of ellipsoid corrections to deflec- 
tions. These contain the values Sa 
and Sf which result when an ellipsoid 
is claanged. These corrections actually 
apply to the geodetic parts of the astro- 
geodetic deflections. 



When it is desired to find a better orientation of the ellipsoid 
as well as to find its size and shape, the observation equations which 
may be written for each point take the following general form; 

(4-1) f = (cp' - cp) + d§ p + df 

= (X^ - X) cos cp + d?^p + d? 7 ^ 

If it is desired only to find a and f, vising data from a meridian 
arc only, the series of observation equations will be as follows: 






1^, 




I 







20 



(4-2) = (cp^' - cp^)+ (d 1^)^ 

^2 = (^ 2 ' “ ^ 2 ^+ ^ c ^2 

***** 

p =(o^-cp)+(d|) 

^n ' ^c‘'n 

2 

Then the least squares solution -inLll involve making = 

minimvun, this giving us the parameters a and f of an ellipsoid which 
will best fit the condition which, in effect, forces the residual de- 
flections to a minimum, as close to zero as possible. Weights can be 
assigned to the observation equations as deemed necessary. 

In general form, the condition can be stated: 

(4-3) I](^^ + V ^) = minimiun, 

2 2 

which is good for any problem. Either or both terms, | and can 

be used according to the problem to be solved. 

In this case, in which astro-geodetic deflections are used, the 
better results are obtained when these deflections have been reduced 
isostatically. 

4.2 Formulas for Displacement and Change of Ellipsoid Corrections. 

There are four well-known methods of computation of the change 
of ellipsoid and displacement corrections to the astro-geodetic 
deflections. These are based on formulas representing different 
derivation approaches to the problem. They are the formulas of Hayford 
( 10 ), the U. S. Array l>lap Service (4), Bomford (2), and Vening Meinesz 
( 33 ) • The derivations of these foimiulas are given in the references 
cited (except (4)), and being beyond the scope of this paper, will not 
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be covered here. 

4.21. Hayford's Formulas. The form of olservation equation given 
in (lO) is as follows: 

(4-4) i = (cp) + 1^ (\) + (A) + n^{^) + o^ (10,000 e^) + (c^'-cp) 



7^ S kg (cp) + 1^ (X) + mg (A) + °2 

+ (\^ - \) cos cp 

in which (cp), (X), (A) are unknown corrections to the initial cp,X, and 

s 2 

A at the origin of the survey, and (]]^) ^<3- (10,000 e ) are unknown 

2 

corrections to the values of a and e for the reference ellipsoid. It 
can be seen that these are similar in form to equations (4-l). 

The change of ellipsoid corrections are: 



(4-5) dg^ = n^ i^) + (10,000 e^) 

= S (ife) + °2 



'c 2 '100' 
100 



N 



in which n^ = - .in ... n 
1 a sin 1" R 



0 cos 






^ -4 



N 

'^l " R 



2 2 

e 0 2 2 

1 + ^ cos cp^ cos Ap 



6(l-e'") 



. 2 
sin cp 



20,000 sin 1" (l-e sin <p^) 



^ g 0 cos Ag 



N 



= “ 

2 a sin 



40,000 sin 1" (1-e ) 
100 



^ (9-cPq)[ 1+2 sin 9^ + 3 cos(cp^+cp)] 



■pr 0 sin Ag 



. 2 
sin cp 



02 = - 



0 sin 






20.000 sin l"(l-e^sin^cp ) 

2 2 1/2 2 ° 

(l-e sin 9^) cos g a, . . 

5-5 — S 8 cos Ap sin A^ 

60.000 a sin 1" (l-e'^) ^ ^ 
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The corrections for displacement (orientation) eire:- 



( 4 - 6 ) (cp) + 1^ (\) + (a) 

" ^2 + Ig (>^) + “2 

q. sin^(X-\ ) 

in which k^ = - 1 + ^ (l + Q) 



= zero 



“l = I ( 1 + Q ) 



2 sin‘ f (Aj - Sj,) 

sin Ag [1 + cos (\ - \^)] 
2 sin^ I (A^ - Ap) 



k^ = - sin cp sin ,(\ - \^) 
Ig = - cos cp 



m^ 



cos cp cos Ag sin (\ - \^) 



sin 






9 = ^ cos= I (fl^ - Ay) 



The purpose of the coefficients kj 1 , ra, n, and o can be stated 
as in this example for one of them: k^ is a numerical coefficient, 

such that if the latitude of the initial point were corrected by the 
amovmt (cp), the change produced in cp* - cp would be k^ (9). 

It can be seen from inspection of the coefficients that for 
points on the same meridian as the initial point, the values are 
k^ = -1, 1^ = 0, and m^ = 0, so that the observation equation would be:- 

( 4 - 7 ) i = (q>) + (3!^) + o^ (10,000 e^) + (9* - 9). 



4.22 Army Map Service Formulas. In an article in Transactions of the 
American Geophysical Union, No. 37 , 195 ^, Chovitz and Fischer 
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reported on the U, S. Army Map Service's determination of the figure of 



the earth from arcs ( 4 ). The procedure is described in some c’etail, 
the essential part regarding correction formulas being reproduced in 
somewhat different form here. The method is designed to be correct for 
meridional arc segments of 1°, but not more than 3° in length, and for 
arcs of parallel of comparable lengbli. 

For meridian arcs the formulas are: 

2 



(4-7) 






8 a , 
a 



- 9^) + 5e2 {l(v9o) -I"" (\ - Vl) 

n=l 



II cos 2cp - III cos 4cp + IV cos 6cp f 
L % ^ra ^raJJ 

in which cp^ = mid-latitude of segment ^ (9^ ^ + 9^) 



CSC 1 " 



(4-8) dEj^ = ~(9o^ “ ” initial orientation of the arc segment. 



In this method the con^jutation proceeds segment by segment, with 
a new initial point for each segment of the arc when the length limits 
given above are reached. Thus the procedure is similar to the one 
described in section 3* 3* This presupposes that all the points in the 
arc furnish data that is adjusted to the same datxun. In contrast, the 
Hayford formulas (and those to be given subsequently) are designed to 
-over a large, but limited area, using only one initial point in the 
computation. The unknowns in formulas (4-7) and (4-8) are 5^, Sa, and 
Be . The Roman numereil coefficients are: 



I 

II 



1 


■7 


17 4 


" ‘ 4 


- 


■ 32 


3 


^ 3 ^2 


^ 45 4 


= 17 


IT % 


5IT ® 




3 


. 9 ^4 
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IV 



3 

SIT 



4 

e 



III 









% 






" Ilk 



• • 








2k 



For arcs of parallel the formulas are: 



) cos cp + <5e^(\ -X ) ■■ - 

'C L a ^ n o' ^ 'no' ^ 



2 . 2 l^■ . 4 X 

(1+e sin cp + e sin cp) 



CSC 1" 



where cp = latitude of the parallel. 

(4-10) d Yo"^^o ~ ^o^ coscp = initial orientation of the arc segment 



2 

The unknowns are >1^, S&, and Se . Solutions for meridians and 
parallels can be done separately or can be combined for a cross-solution 

(4). 

4 . 23 . Bomford's Formulas. The formulas for change of ellipsoid 
given in Bomford (2) p. 127-130^ actually contain corrections for both 
displacement and change of ellipsoid. They are developed here in 
separated form. These formulas are used in Bomford's version of the 
area method^ which supposes a "homogeneously computed national survey" 
covering a 20° by 20° area containing 50 astro-geodetic stations. 

Bomford states that it will generally be best to ignore values of 
deduced from azimuth observations (these were omitted in the description 
of Hayford's method above), and to use only those based on longitude 

(2) p. 365 . 

The change of ellipsoid formulas as given by Bomford (2) p. 129, 
365 , are as follows: 

(4-11) d| =- ^ 'l/a)sir. u cos CD (U sin u^ + \ cos u^) 

+ sin u sin ~ (l/a)cos u(V sin u^- U cos u^) 

CSC 1" 



+ sin 2u Sf 






-<r. 
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(14--12) d.y( = - ^ (l/a) sin a> (U sin + V cos u^) - cos to esc 1" 

where U = -aS| - a<5f sin 2 u 

^o o 

r r 2 

V = W “5a + aof sin u 
o o 

u = reduced latitude; tan u = (l - f) tan cp 

CD = longitude measured east from origin 

= height of origin in corrected ellipsoid 

above that in reference ellipsoid 

Si and are in radians. 

^o <o 

Substituting the values of U and V in (4-ll) and (4-12), and 
separating the expressions into two parts: d| = d|^ + d|j^ 



d>7 =d>^j^ + d>^jj, we obtain: 



(4-13) dl^, = 



(sin u cos CD cos u - cos u sin u ) 5 a/a 

' o o' ' 



^si] 



.. 2 



+ sin u cos CO sin u cos u + cos u sin u 

o o o 



(l + cos^u^) - sin 2uj Sf 



CSC 1” 



(4-l4) df^ = (sin u cos cd sin u + cos u cos u ) 56 " 

^ ' o o'^ ^*o 

- sin u sin cd SW'- (sin u cos cd cos u 

•o' o 

- cos u sin u )(N /a) esc 1" 

o' o' ' 

(4-15) <^70 = l^sin CD cos u^( Se./e) + sin co sin^u^ cos u^ 5fj esc 1" 
(4-l6) d>^jj = sin CD sin u^ 5||^ + cos cd S>^” 

- sin CD cos u (N /a) esc 1" 

o o' ' 

4.24, Vening Heines z' Formulas. The formulas here referred 
to are those published in 1950 ( 33 ) } sind which were developed for the 
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purposes of (a) describing the relationship between the variation of 
the vertical deflection components at the origin of a survey and the 
changes resulting in the deflections at any arbitrary point in the 
survey, the variation having been caused by a parallel displacement 
of the reference ellipsoid with respect to the geoid; and (b) 
describing the changes in the deflections of the vertical brought 
about by a change of reference ellipsoid. 

The complete formulas referred to in (a) above, are as 

follow: - 

(4-17) d|j^ = [cos(cp-cp^)-2 sin cp sin sin'" i(\-\^)] 
r 

+ " sin cp sin 

1 2 1 
+ ^ [sin (cp-cp ) - 2 sin cp cos cp sin ^ (\-\ )]&N 
R ' o' ^o 2 . o' o 

R r 

(4-l8) dT)_ sin cp sin ( )6| + — cos (X-\ ) m 

' ' 'D r ^o ' o' ^o r ' o' o 

- i cos cp sin ( X-\ ) 6N 

r ^o ' o ' o 

where R , R = radius of curvature in meridian at initial point P 
and at point P, respectively 

r r = radius of curvature in prime vertical at initial 
o, 

point P^ and at point P, respectively (r is used 
here, rather than the usual N, which is reserved 
for geoid distance). 

The above form is that given by Heiskanen in (13) p. 302, and 

is better suited to the use of tables than the form in the original 

2 3 

paper (33), which uses a(l-e )/w^ rather than R, and a/W rather than 
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2 2 ^ 

r, where W = (l-e sin cp) . 

Vening Meinesz points out that his formulas are similar to 
those developed by Helmert (37 )> but include 5N^ terms as essential 
differences , 

If the line P^P is long, the complete formulas above must be 
used. However, for nets of small east -west extent, terms 
proportional to the squares of can be neglectedj and for 

nets of small extent in latitude, the ratios 

r^/r are very nearly 1 (within 0.01), and the formulas (4-17) and 
(4-l8) become: - 

(4-19) d|_ = cos (cp-cp )6| + sin cp sin (\-\ )5t) 

' D ^o ^ ' o o 

5N 

+ sin (cp-cp^) ^ 

(4-20) dT)j^ = - sin cp^ sin (^.-X^) 5|^ + 6 t)^ 

5No 

- cos cp sin (X-\ ) 

According to (13), the formulas (4-19) and (4-20) apply when 

P P< 300 km. in east-west direction and < 1500 km. in north -south 
o 

direction. However, Vening Meinesz (33) states that if it is 
desired not to neglect quantities greater than 1/300 of the 
deflection of the vertical, these formulas can be used for nets 
extending up to 4° 4o’ in east-west direction euid to more than 20° in 
latitude, from the origin P^. Furthermore, the complete formulas 
will seldom require that the ratios between R, R^, r, and r^ be any- 
thing but 1 even for large nets, to attain the required accuracy. 






r 










= - m € 
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The complete formulas referred to in (b) above, for change of 
ellipsoid, are here transcribed; - 

(4-21) d|^ = |[sin (cp-cp^) - 2 cos sin qj sin^ (\-X^)] 5P 

1 X 

- 4 cos cp cos ^ (cp+cp^) sin ^ (cp-cp^) 5f 

2 

+ [(2-3 sin q) - cos (cp+cp^)) sin (cp-cp^) 

^2 2 2 X % 

- 2 cos q)^ sin q (2+sin - 3 sin q)sin ^ 

- sin 2q^ sin - (q-q^) - 2 cos q cos g- (q+q^) 
sin I (q-q^)] f&p + [ | sin 4 q^ (sin^ | (q-q^) 

+ sin q^ sin q sin ^ (X-X.^)) - sin q^ sin ^ (q+q^) 

3 XX 

sin 2 (9-Tq) - sin q^ sin g- (q+q^) sin ^ (t-9q) 

- 4 cos q cos I (q+q^) sin | (q-q^) 

+ 4 sin 2 q cos^ ^ (q+^p) sin ^ ^ (q-q^)]f6fj cscl" 

(4-22) dT)^ = |- cos q^ sin (X-X^) 5p 

+ cos q^ sin (X-X^) [sin (q-q^) sin (q+q^) f&P 

+ ^ tan q^ sin 4 q^ f6f]^ cscl" 

Sa 2 

where 5B = — + sin q 6 f 

The formula (4-21) is quite complicated, but can be 
simplified for systems not of world-wide extent. When the terms f6f 
and fSp are assumed to have a magnitude of 1;6,000,000 or less, the 
terms containing the squares and products of (q-q^) and (X-X^) may be 
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neglected in the factors by which the terms are multiplied. When 
these angles are such that the distance from the origin is 600 km. , 
the sqxiares and products amount to an order of magnitude of 
1;600,000,000 (33). 

The simplified formulas ares- 

(I+-23) d|^ = |^[sin (cp-cp^) - 2 cos cp^ sin cp sin^ | (X-X^)]6 p 

- k cos cp cos ^ (cptcp^) sin ^ (<P-<Pq )5 f 

- (2 + ^ tan cp^ sin kcp^) sin (9 -9^)f6f| cscl" 

( 4 - 24 ) dn = [- cos 9 sin (X-\ )6 p 
' ' c o o 

+ ^ sin 9^ sin 4 9^ sin (X-X^) f6f] cscl" 

4.3 The Use of Gravimetric Deflections of the Vertical. 

Probably the most accurate method of determining the 
equatorial radius from deflections of the vertical, involves a final 
correction by application of the gravimetrically determined 
deflections of the vertical ("absolute" deflections). 

The following analysis could have been given at the 
beginning of this chapter, but it was felt better to place it here, 
thus keeping the complete picture intact. 

Rice (25) p. 2, states that "a fundamental geodetic problem is 
to arrive at a reference ellipsoid of such dimensions and having such 
orientation that the deflections of the vertical are as small as 
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possible in the aggregate." 

An ideal reference ellipsoid would be centered at the center 
of mass of the earth and its polar axis would coincide with the mean 
axis of rotation of the earth. This is what Veis (32) p. l6, calls 
the Terrestrial Ellipsoid. A further requirement which would be 
placed on an ideal reference ellipsoid would be that its surface best 
fit the surface of the geoid. 

The origin or initial point of a geodetic system is 

a 

generally assigned the astronomic latitude and longitude, Cf/ and 
observed at the place. The geoid distance is generally considered 
zero and no errors in azimuth are assumed. The effect of making 
(p = <p*, £ind \ = \*at the initial point is a parallel displacement of 
the reference ellipsoid’s system of coordinates with respect to those 
of the terrestrial ellipsoid and the geoid. At another point, P^, in 
the same geodetic system, which is displaced from the geoid, the 
normals do not coincide, the result not only of the displacement, but 
of the wrong computation surface, which is not the same size or shape 
as the true earth spheroid. The resulting meridian angle between cp*' 
and cp is the astro-geodetic deflection (cp’-cp) in the direction of the 
meridian. (This discussion is generally in terms of this meridian 
component, for convenience, though similar remarks can be made 
concerning the prime vertical component, as well as the total 
deflection). 

' The meridian component and the prime vertical component of the 
astro-geodetic deflection each consists of three angles in the same 





4 




9 
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plane, with respect to the normal to the mean or smoothed geoid. If 
we can correct this deflection in such a manner that both components 
are reduced to zero, what would remain would be the normal to the 
mean geoid. In doing so we have to change the size, shape, and 
position of the reference or computation ellipsoid so that our point 
coincides (as nearly as possible) with the corresponding point on 
the mean geoid, and so the new ellipsoid is in the position of a 
terrestrial ellipsoid. The three eingles ares- 

(1) Correction for displacement, 

( 2 ) Correction for change of ellipsoid, 

(3) The absolute deflection of the vertical, or the angle 

between the true vertical and the normal to the 
mean geoid, or if preferred, the terrestrial 
ellipsoid of best fit. Since this is obtained 
gravimetrlcally, we call its two components | 

S 

and . 
g 

See figures 3 and 4 for pictorial descriptions of this. 



32 

AST/ZO- G£OZ>ET/0 




Fig. 3 

MEAsf 

T/ec/£. i/E/ZTttTAl. \/£/ZT/CAL. 




Pig. h 



33 



It should be noted that when the gravimetric deflections at 
the initial point are known, two of the unknowns of the previous 
problem are now known, namely 6|^, and These are equivalent to 

the following 



(i^-25) 



where 



5 So = - (ta - Sg)o = 

5 1o ' - (la - \>o ' (ig - ’'a'o 

(Sa>o = ('»o - 

(''a>o ■ (^o • ‘Pq 



The gravimetrically computed quantities | and q , are in- 

S S 

dependent of the size of the reference ellipsoid (13) p. 239> even 
though they depend on gravity anomalies derived from the gravity 
formula employed. 

In the method (described in section 4.1) of determining the 
dimensions of the earth ellipsoid by means of correcting the astro- 

2 

geodetic deflections with d|j^ and d|^, and subsequently reducing 
to a minimum, we have one effective method of solving the problem. 

However, we can go one step farther, as indicated in the preceding 

paragraphs, and make the final correction by subtracting the value | 

g 

from the astro-geodetic deflections which have already been corrected 
by d|jj and d|^. 

Thus our problem will be solved by a least squares solution in 



which 
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2 (| - I )^ = minimum 

g 

p 

2 (ti - T] ) = minimum 

g 

In this way we are much closer to the mean geoid at the start. 

The result shovild he a terrestrial ellipsoid which best fits the 

surface of the geoid, providing, of course, that we have sufficient 

points from a large enough area for our computation. Theoretically 

the problem can be solved when astro -geode tic and gravimetric 

deflections are known at just two points, one at each end of a long 

•% 

arc, preferably as long as the earth^s radius. Naturally, the more 
points at which this data is known the better the solution. Heiskanen 
has suggested this method in a nviraber of publications; see par- 
ticularly ( 12 ), and ( 13 ) p. 278 - 279 * 



Chapter 5* Flattening from Gravity Anomalies. 



A value for flattening can be gotten from the solution by the 
astro -gravimetric method described in section but if this is done 

the gravity anomalies used in computation of the gravimetric deflect- 
ions are no longer correct since these anomalies are based on the 
adopted gravity formula, which in turn is based on a particular 
ellipsoid with its own flattening. This gives a clue as to how the 
flattening may be corrected fi;om the anomalies themselves. See (13)^ 

p. 279 . 

According to the theorem of Clairaut (13) p. 52, the 
following is true (extended form):- 



(5-1) 





f 



17 
• I? 



mf 



where 



7 f 7y are normal gravity at pole, equator, 
P 

m = (centrifxigal force at equator )/7' 

Jbj 

f = flattening 



resp. 

2 

CD a 



The basic form of the gravity formula giving normal gravity 7 
for any latitude q) on the ellipsoid is (neglecting the longitude term) 



(5-2) 7 = (1 + p sin^cp - 6 sin^ 2 cp) 

5 12 

where € = ^ mf - ^ f 

The problem is to correct the parameters 7 , p, and e, to the 

Jlj 

gravity formula where these parameters have been computed for a 
particular reference ellipsoid. Thus if p and e are corrected, a 
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correction to f is automatically found. The correction to e (based on 
currently used ellipsoids ) is always so small, we may neglect its 
computation. 

Using, for example, the International Ellipsoid, the following 
values are used as a basis 



7 978.049 gal 

p = 0.0052884 



a = 6 378 388 m. 



f = 1/297 = 0.0033670 
m = 1/288.36 = 0.0034678 



6 = 0.0000059 



For obtaining a corrected value of 7, we have;- 



where 




in which x' = unknown correction to 7 , in mgals 



y* = unknown correction to p, in mgals 



Before correction, the gravity anomalies, based on the 



ellipsoid used and isostatically reduced, are:- 



^ - 7 



After correction they are:- 



V = Sq - 7" = Sq “ (7 + A7) 
(gQ-r) - A7 = - A7 



Then, we can write error equations for each point where 



gravity observations have been mades- 




If we have enough such stations, spread out over a large 
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enough area, we can solve the prohlem effectively by the method of 
least squares in such a manner that;- 

2 

Z V = minimum. 

An example of this type of solution made by Heiskanen, in- 
volved using anomalies from over 1500 squares of 1° by 1° dimensions 
( 13 ) P» 77 * The resulting parameters were:- 

= 978,0451 

p =0, 0053026 

€ = 0.0000059 

The flattening corresponding to the corrected P was:- 



f = 1/298.2 



•/m 



j 



Chapter 6. Gravimetric Deflections of the Vertical 
6,1 General. 

The gravimetric method, of which the subject of Chapter 5 is a 
part, also produces what are sometimes referred to as "absolute" 
deflections of the vertical, obtained from gravity anomalies. The 
absolute or gravimetric deflection of the vertical at a point is 
defined as the angle between the normal to the geoid and the normal to 
the earth spheroid at the point. — Veis (32) uses the term "mean 

•a 

geoid" to describe the surface of the earth spheroid in the area of a 
point, meaning a smooth surface striking a "mean" between undulations 
of the geoid, — The discussion in Chapter 4 (section 4.3), describes 
the use of these deflections in determining the size of the earth 
ellipsoid. This chapter will present a brief outline of the method of 
determining the deflections from gravity anomalies. 

To begin with, the gravity anomalies /2g, are caused by dis- 
turbing mass anomalies, and with this in mind and with the aid of 
fig. we can easily see how local absolute deflections occur. 




Fig. 5 
38 
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The mass surplus tends to "pull" the plvunb line toward it, and 
conversely the mass deficiency tends to repel it, thus causing the 
deflection S = "P* - The sign of this deflection is apparent from 
this expression which shows that | is always referred to the 
astronomic latitude, cp', of the point in question. The sign conven- 
tion for gravimetric deflections is the same as for astro-geodetic 
deflections. Referring to figure we wee that when cp > cp', | = 

(cp' - cp) < 0, and when cp < cp' , | = (cp' - cp) > 0. In U.S. practice \ 
is considered positive westward# and in European practice positive 
eastward. In either case, when \ > \' , t) = (\' - x) cos cp < 0, and 
when X < X' , f = (X' - X) cos cp > 0. 

From a theorem developed by Stokes, Vening Meinesz formulated 

expressions for absolute deflections of the vertical in 1928. Stokes' 

Formula expressed the geoid distances N, as follows 

2tt w 
R 



(6-1) N = 



2 7T G 






f(t) sin t ^ d t 



in which \|c = angular distance of a circular ring of width d\lr. See 
fig. 6. 

Ag = average gravity anomaly in ring, having the limits t and 

t + dt. 



R = Mean radius of the earth. 

G = Mean gravity. 
f(t) = I S(i)o 

S(\[r) = CSC i t + 1 - 6 sin ^ \|r - 5 cos t 



(6-2) 



4o 



1 2 1 

- 3 cos \lr 1 n (sin ^ 'if + sin ^ t) 

which is the original Stokes Function, 
sin t d\lr dA = surface area element for unit sphere. 



POLE 




Vening Meinesz got his formulas by differentiating the Stokes 
formula, as indicated in fig. 7> in the form:- 

N = 2 - ?Vr /^■(^)^dS 

s 




Upon differentiating, the following are obtained:- 







•• 




— • 










I 



4l 



(6-4) 


d|" = 


5 N 


CSC 1" / 


f(t) 


cos A z:^g d S 


" S X 


2 TT G R^ 

S 


(8-5) 


dT|" = 


5 N _ 


CSC 1" j 


f(t) 


sin A d S 


" 5 y ' 


2 IT G R^ ^ 



2 

where dS = R sin i|f d ijr dA, and 

cos ^ 

(6-6)f (ijr) = = - |-x- + 8 sin ijr - 6 cos |i|f 

^ 2 sin ^ 

- 3(1 - sin CSC ijf + 3 sin In [sin ^ (1+ sin ^)] 

Application of these formtilas to practical use requires 
further manipulation. They are used in different form, whether one is 
going to conqpute the effect of areas beyond « 3°, or within the area 
circumscribed by a circle of radius \|f 3° (3° is an arbitrary dividing 

line, which, though close to the limit, has been chosen for convenience 
(31) p* 98 )* These different forms are described as the "Circle-Ring" 
Method, for use up to ijf 3°, and the "Squares" Method, for use 
beyond. Still another form is used for the small inner circle close to 
the computation point, having a radius r^, generally, of 5 km. or less. 
Figure 8 shows this division of methods pictorially. For further 
references, see (13) p. 63-70, 257-264, (15)« 

6.2 The Squares Method. 

This method, used for the earth represented as a spherical 
surface, and for the area beyond a radius of about 3° from the comput- 
ation point, makes use of mean gravity anomaly values for each 1° by 
1° square (spherical trapezoid) between « 3° and \|r ~ 20°, and for 
each 5° by 5° square beyond 20° aud around to the antipode (31)* 


















# ^ i I # ^ 4 #• f - 




These should be either free -air or isostatic anomalies, the latter 

being best. If Bouguer anomalies are used they must be corrected to 

obtain free -air anomalies. To adapt formula (6-4) to the squares 

2 

method, we first replace dS by R sin i|r d ilr d A, obtaining; - 

(6-7) d|" = / f (t) cos A Ag sin i|r d t d A 

s-' 

Then using the sqmre element dq on the unit sphere, R = 1, 
(l4), we obtain; - 

(6-8) d|" = 2 ^" - ^- J"f ' (ilf ) cos A /:;« d q 

The expression for dii" being the same with sin A replacing 

cos A. 

2 

When ^ is in milligals, and when G is taken as 9T9»^ cm./ sec. 
the term before the integral sign will be; v 



y 



h3 



c 



ir 



CSC 1^’ 

2 7T G 



= 0'.'0335 



One procedure is to perform numerical integration summing up 
the effect of each square over the earth as indicated in the following 
formula; - 



its mean latitude. The subscript q, where occurring, signifies mean 
values of the quantity for the square in question. 



electronic computing machine, but manual graphical solution, though 
more time consuming, is made easier by use of maps on polar stereo- 
graphic projection. With this projection templates can be easily 
prepared for use on the large area, small scale maps required, since a 
feature of the projection is that all circles on the sphere project as 
circles on the map, so azimuth lines (which are great circles) and the 
circular rings (small circles) can be easily drawn. Necessary 
formulas for preparing such templates may be found in (13) p. 259-260, 
and in (31) p. 98-100* 



been computed by Sollins at close intervals for different angular 
distances t, from the computation point. These tables have been 
published by Sollins (29) and are very useful in the solution of the 
problem. 



(6-9) d|" = c"Zf'(t) cos A ^ q 



where q is the area of a 1° by 1° or 5° by 5° square and depends upon 



The squares method adapts readily to solution by high speed 



The functions f'(\|r), f*('lf) sin and F'(\lr) sin t d \]f, have 




kk 



6o 3 The Circle -Ring ^fethod. 



In the relatively close area around the computation point 
(within \lf ~ 3°)^ where a detailed gravity anomaly map is available, 
formula (6-7) is used in the following form:- 



( 6 - 10 ) 

d|" 



CSC 1" 
2 7T G 




\lr+d\|f 




sin d 



CSC 1" 
2 TT G 



(sin - sin 



\jf+d\lr 




f*(t) sin \lf 



d ^ 



For d T]", substitute (cos - cos Ag) for (sin A^ - sin A^). 

Formula (6-10) gives the effect of a circle-ring compartment bounded 

by ^ and \|r + dijr and by A^ and A^, and having a mean anomaly on 

the total component | 

r ® 

Values of ^ (\|r ) sin \Jf d are tabulated in Sollins Tables. 

4=10m \lf4d\lf 

Thus for the value of a particular y' f* {'^) sin \lr d \lr, the tables 

would be used as indicated in the following expression: - 
t+d\lr 






(\(r) sin \lf d t 



\lr=10m 



-I 

T|r=10m 



f' (\(r) sin \lr d \lr. 



Figure 9 gives pictorial explemation of the components of the 
method and corresponding templates. 

Any system of azimuth sector and ring zone intervals can be 
used in this method, but a system which yields equal effect for each 
compartment would obviate much unnecessary work. Kasansky developed 
a system in which the sine differences are l/8 (32 compartments), and 



J\/OR.TH 
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i + 
A 



P = computation point 
\lf = inner radius of compartment 
dt = outer radius of compartment 
Ag = bounding azimuth lines 

= mean gravity anomaly in mgals, read from 
iso -anomaly lines on map 
Fig. 9. Components of Template for Circle- 



Ring Method 
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the ring radii, beginning with r^ = 5«0 km., follow the geometrical 

progression r = 1.2T0°r . Farther from the station, the ratio is 
no 

slightly smaller than 1.270 because of the earth ^s curvature. Tables 
of these Kasansky Circle -Rings may be found in (13) p. 267, and in 
(16) p. 17. The effect of each compartment is O'.'OOl radial deflection 
effect per 1 mgal 

D. A. Rice, of U.S.C.& G.S. , published a system of Circle- 
Rings (25), in which each ring was divided into 36 equal azimuth 
sectors of 10° aperture. Her§ also, the radial effect for each com- 
partment for 1 mgal ^ is O'.'OOl. Rice^s method for computation of the 
radii involved using the formula for the deflection effect along any 
radial lines- 



(6-11) 


<30" = 2% ' g ~' J ^'(^) sin d t ^ d A 

s 


SO 


d|" = d 0" cos A 




dV' = d 0" sin A, and 


(6.12) 


CSC 1" ^1 



t-, 

where I. is obtained from Sollins^ Tables (29)* 
^o 



Then for AA = 10° = 0.17452 radians, and for G = 98l>000 mgals, 
formula (6-11) becomes :- 

^1 

d0" = 0.0058405484 I 

^o 

For d0" = O'.'OOl and z^ig = 1 mgal. 




= 0.17121680 






t 



* « ^ • 1 1 .— * 
» t • 4| - t ™ 









I 




i«ii»» * •‘Jfc •' 
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i \]f 

Alg = = 0.17121680, etc. 

Then selecting interpolated values of \|r from Sollins* table 
giving Al, the circle-ring table was prepared. Table I is Rice's 
table reproduced. 



Table I 



Rice Template Radii for Gravimetric Deflection of the Vertical 



For radial deflection effect of O'.'OOl -and mean com- 
partment anomaly of 1 mgal; angular aperture 10 . 





(n 


= zone 


nvimberj r = 


inner 


radius in 


km. ) 




n 


r 


n 


r 


n 


r 


n 


r 


1 


0.119 


14 


1.099 


27 


10.15 


39 


77.97 


2 


0.l4l 


15 


1.304 


28 


12.05 


4o 


92.22 


3 


0.167 


16 


1.547 


29 


l'!-.29 


4l 


109.0 


k 


0.198 


17 


1.836 


30 


16.94 


42 


128.7 


5 


0.235 


18 


2.179 


31 


20.09 


43 


151.9 


6 


0.279 


19 


2.5^ 


32 


23.83 


44 


179.1 


7 


0.331 


20 


3.068 


33 


28.25 


45 


210.9 


8 


0.393 


21 


3.641 


34 


33.48 


46 


248.0 


9 


0.467 


22 


4.320 


35 


39.67 


47 


291.2 


10 


0.554 


23 


5.125 


36 


47.00 


48 


341.2 


11 


0.657 


24 


6.081 


37 


55.66 


49 


399.0 


12 


0.780 


25 


7.216 


38 


65.90 


50 


465.5 


13 


0.926 


26 


8.560 












With the 


use of 


circle rings 


such 


as Rice's, 


the problem of 



computing the gravimetric deflections becomes simple, though laborious. 
All that is necessary is to estimate the mean anomaly from the contours 
in each compartment, make svimmation of for each azimuth sector, mul- 
tiply by 0.001 cos A or 0.001 sin A , and the effect for that sector is 

m m 

solved. Then the effects of all 36 sectors are added together for total 
effect. This may be further simplified when one takes note of 
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the fact that mean azimuths in 150° opposite sectors will have the 
same cosine or sine, but with opposite sign. One must be careful with 
signs, since both the azimuth of the sector and the anomaly gradient 
contribute to the sign. Additional information on this aspect of the 
problem will be given in the explanation accompanying the computations 
made in Chapter 7* 

In the event gravity material is poor, which is too often the 
case (and one has to be a good judge concerning when poor material is 
to be considered good), Heisk^en (13) p. 268, recommends combining 
two zones and two sectors into one four-compartment compartment, the 
effect of which then would be O'.'OOU per ragal. 

6.4 The Inner Circle. 



The circle -ring method, although it can be applied beyond 
= 3°, should not be used very close to the computation point since 
the term esc ^ i|r in the value of f'(t) sin t, in formula (6-10) 
approaches infinity. Rice (25), suggests a lower limit of 100 meters. 
However this is usually too small because of the sparsity of gravity 
anomaly values, except in isolated instances, within any given 100 
meter circle. 

The inner circle, being small, is considered to have a 
uniform anomaly gradient, so the anomalies read from the cardinal 
points are used in the following expressions, the derivation of which 
may be found in (31) p. 93, 100-101:- 

(6.13) dS" . X = 0V105 
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where 



r = radius of inner circle in km. 



o 



6 = Zigg - ZNgjj in mgals (N. hemis. ) 



5 (^)y = in mgals (European practice) 



&x 




by = 2 r 



o 



Substituting these where appropriate, formulas (6-13) become; - 

(6-l4) d|" = 0V0525 (^g ^j^) 

dT)" = 0V0525 (^ - ^g) 

See figure 10 for further clarification of this "one- 
gradient" method of computing the effect of the inner circle. 



Rice used a three -gradient method rather than one -gradient as 
above, in computations for l6 points in the U.S. (26). Taking 






•^Sw 










Fig. 10. Inner Circle - One-Gradient Method 



anomalies at NE, SE, SW, and NW points, as well as at the cardinal 
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points of the circle, three N-S gradients and three E-W gradients were 
used, with weights of one -half given to the outer gradients and a 
weight of one to the center gradient. Figure 11 illustrates this. 

This method should he used only when detailed gravity information is 
available inside the inner circle. 







Sx rs-A/) 

<Sy (tV-£) 

Sx CSw-A/W) 

(SE-NE) 
Stj (S\Af-SE) 
(aj\a^-aie) 



2. r. 



/. 4-14. 



Fig, 11. Inner Circle - Three -Gradient Method 



^ ^ 

d^"(S-W) = 0V105 — ^ r Weight = 1 

o ^ o 



d|"(SE-ME) = 0V105 ■ ■■ Q - x ■ ' Weight = | 



/Vg .. /\g 

d|"(SW-nw) = 0V105 — ^ ^ r Weight = | 

OX O' ^ 

£ Weights = 2 



Combining these, using weights shown and values of 6x given 
in fig. 11, we obtain the following expression, given in (31) p. 101:- 



(6-15) d^-' = 0,0262 (Z^g- dgjj) + 0.0186 
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This is for the northern hemisphere. For the southern 
hemisphere, Just reverse S and N where occurring. The other component 
of the deflection is similars- 

(6-16) dti” = 0.0262 + 0.0186 

This is for European practice, in which X and are considered 
positive eastward. Where the practice is opposite, W and E should he 
reversed where occurring. 

As may be deduced from the opening paragraph, the smaller the 
inner circle, down to certain limits, the more concentrated the gravity 
information must be, Vening Meinesz (13) p. 70, states that the error 
incurred using a circle with r^ = 30 km. is less than 1^. A good 
round number seems to be 5 km. , and for good results there should be 
from four to eight points of known ^ within the circle (13) P* 267. 



Chaptei f. Computations 



7.1 Dimensions from the 98 th Meridian in the United States. 

(From Astro -Geodetic Deflections). 

Using data from U.S.C. & G.S. Special Publication No. 229> (7)> 
consisting of astronomic and geodetic latitudes eind meridian com- 
ponents of the astro-geodetic deflections of the vertical for 30 
o o 

points between 97 W and 99 W longitude, extending from North Dakota 
to Texas, values for a and f will be determined. The astro-geodetic 
deflections as published in the publication named have been reduced 
according to the Hayford isostatic assvimption for a depth of com- 
pensation of 113*7 km. 

2 

The problem to be solved is E | = minimum, according to the 

theory of least squares. The Vening ^feinesz displacement and change of 
ellipsoid formulas will be used. To adapt them to practical use, terms 
containing the same \anknowns must be collected and combined. 

Prior to solution, a test problem was solved to determine what 
terms may be disregarded and whether or not the long or short formulas 
should be used. The maximum difference in longitude is about 1°, and 
terms containing longitude were determined to make negligible con- 
tribution to the residual deflection. Also, for the arcs of limited 
length (Acp~13V3^ in this case) encountered in this problem, the 
difference between residuals computed by long and short formulas, 
turns out to be negligible. Terms containing 6N will also be 
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neglected to simplify the problem. 

The Vecing Meinesz formulas, given as (4-19) and (4-23) in 
section 4.24, are here repeated with terms omitted as noted above. 

Displacements- d^j^ = cos (cp^ - cp^) 6^^ 

Change of Ellipsoids- d^^ = p^6P - PgBf - p^f S f 

8a 2 

where 8P = — + sin cp 5 f 
a ^o 

Substituting the foregoing expression for SP in d|^ and collecting 
terms, we obtains- 

5a 2 

d^c = Pi ^ Pi " ^3^^^ 

Pi 2 

= __ 5a + [p^ sin <Pq - Pg - P3P] Sf 
in which p, = sin (cp - <p ) 

I'l \t-jj -r^y 

p_ = 4 cos qp cos ^(cp +cp ) sin ^(<p - cp ) 
^2 ^n 2'^n ^o'^ 2 n ^o"^ 

O 

P3 = (2 + ij: tancp^sin 4cp^) sin (cp^ - cp^) 



The observation equations will take the form given below, with sub- 
script 0 signifying the initial point, Meades Ranch, Kansas, and sub- 
script n signifying any other point in the arc. 



5 = d^„ + di +(<+ '- cp ) 

n D c n ^n'^ 

For convenience we transform the unknowns and the coefficients 
as indicated heres- 



8 a 

In = (h)B|^ (j) 1000 + (k) 10,000 8f + (l) 



= hx + Jy + kz + 1 
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where x = 6| 



and 



y - 



z = 



o 

Sa 



1000 

10,000 6f 



unknowns . 



cos (cp - CD ) 

'^n ^o' 



h 

j = sin (cp^ - cp ) 



1000 CSC 1" 



k = 



[sin (cpn-9Q)[ sin^ cp^ - (2-t^ tan cp^ sin 4 cp^)f] 
- 4 cos.cp_^ cos I (<p„t<p^) sin | (<P„-<P„)] 







• I JP W f 

1 

*8 - 



n 



Table II - Basic Data - 9Q~th Meridian 



Point 




Latitude Astro-Geodetic 


n 


Location 






Deflection 








n 


(cp» - cp) 


496 


W.Dak. 


46 ' 


^ 02 ' 17" 


'N. 1'.'17 


571 


S.Dak. 


45 


28 16 


1.06 


566 


t» 


44 


53 53 


0.78 


569 


M 


44 


00 52 


-2.93 


570 


II 


43 


43 24 


-1.43 


567 


II 


43 


18 48 


0.96 


412 


Webr. 


42 


25 29 


3.88 


407 


II 


41 


35 44 


-1.29 


413 


II 


4 o 


46 35 


-2.93 


235 


Kans . 


39 


46 35 


-2.03 


Meades R. 


II 


39 


13 27 


0.00 


238 


II 


38 


43 47 


-1.29 


245 


II 


37 


51 38 


-2.51 


24 o 


II 


37 


14 28 


- 0.40 


522 


Okla. 


36 


4 o 08 


0.20 


524 


II 


35 


56 49 


1.88 


516 


II 


35 


16 27 


1.71 


521 


II 


34 


38 53 


1.98 


602 


Texas 


33 


37 21 


-1.28 


620 


II 


32 


48 10 


-0.91 


639 


II 


32 


15 


6.78 


605 


II 


31 


39 26 


2.20 


59 ^^ 


II 


30 


59 ^ 


7.65 


592 


II 


30 


16 21 


-1.11 


635 


•1 


29 


42 49 


- 4.05 


652 


II 


29 


22 16 


2.14 


627 


II 


28 


52 38 


-0.59 


634 


II 


28 


18 06 


-0.68 


588 


II 


27 


44 35 


1.44 


630 


II 


27 


4 o 16 


4.91 


darkens 


1866 Ellipsoid:- 




a = 


8 , 378,206 m. 








f = 


1/294.98 










0. 00339006 



Initial Points- Meades Ranch 

cp^ = 39°13'26V686 n. 
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Table III - Observation Equations - ^Qth Meridiaxi 



Point 


h 




k 




n 


(x) 


(y) 


(z) 


= (<P' ”9 



J+96 


0.993 


3.837 


-1.544 


IV 17 


571 


.994 


3.519 


-1.450 


1.06 


566 


.995 


3.197 


- 1.348 


0.78 


569 


.997 


2.701 


-1.178 


-2.93 


570 


.997 


2.535 


-1.118 


-1.43 


567 


.997 


2.306 


- 1.033 


0.96 


kl 2 


.998 . 


1.806 


-0.835 


3.88 


ko 7 


.999 


1.338 


- .637 


-1.29 


413 


1.000 


0.876 


- .429 


-2.93 


235 


1.000 


.312 


- .158 


-2.03 


Meades Ranch 


1.000 




-- 


0.00 


238 


1.000 


- .279 


.146 


-1.29 


245 


1.000 


- .770 


.414 


-2.51 


24 o 


.999 


-1.119 


.613 


-o. 4 o 


522 


.999 


- 1.442 


.803 


0.20 


524 


.998 


- 1.849 


1.051 


1.88 


516 


.998 


-2.228 


1.290 


1.71 


521 


.997 


-2.580 


1.520 


1.98 


602 


.995 


-3.157 


1.911 


-1.28 


620 


.994 


-3.617 


2.367 


-0.91 


639 


.993 


-3.920 


2.457 


6.78 


605 


.991 


-4.258 


2.710 


2.20 


594 


.990 


-4.628 


2.993 


7.65 


592 


.988 


-5.032 


3.310 


- 1.11 


635 


.986 


-5.353 


3.567 


- 4.05 


652 


.985 


-5.534 


3.716 


2.14 


627 


.984 


-5.808 


3.944 


-0.59 


634 


.982 


-6.128 


4.214 


- 0.68 


588 


.980 


-6.437 


4.481 


1.44 


630 


.980 


-6.477 


4.516 


4.91 



i 
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Table IV - Least Squares Solution - ^th Meridian 



Forward Solution: - 










X 


y 


Z 


1 


Z 


X +29.62 

-1 


-^'■ 7.^3 

+1.6oi2d3 


+ 35-76 

-1.207292 


+15 08 
-0.509115 


+33.0? 

-1.115125 


y 


+398-56 
- 75-95 


-250.22 
+ 57.26 


-88.55 

+24.15 


+12.36 

+52.89 


( 1 ) 


+322.61 

-1 


-192.96 - 6 !+. 4 o 

+0.598122 +0.199622 


+65.25 

-0.202257 


z 




+160 . 48 

- 43.17 
-115.41 


+58.88 

-18.21 

-38.52 


+ 4.90 
-39.88 
+39.03 


( ) 




+ 1.90 
- 1 


+ 2.15 
-1.131579 


+ 4.05 
-2.131579 


1 






+223.62 

- 7.68 

- 12.86 

- 2.43 


+209.03 

- 16.82 
+ 13.03 

- 4.58 


[ll] 






+200.65 


+ 200.66 


Back Solution: - 










z 


Check 


y 


Check 


X Check 


1 -1.131579 

z 

y 


-2.131579 +0.199622 
-0.676822 


-0.202257 

-1.274944 


-0.509115 -1.115125 
+1.3661146 +2.573438 
- 0.764132 -2.365417 


z -1.131579 

Z(check-U) 


-2.131579 

- 1.000000 


-0. 477200 


-I.4772OD 

- 1.000001 


+0.092899 -0.907104 
-1.000003 



Standard Error of single observation: - 

= ± J ^ ^ 3 ol 3 ^ = ± nT 7.^^3145 =i 2'.'726 
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Table V - Computation of Correlate Numbers - 98th Meridian 





l(x) 


I(y) 




l(z) 


(x) 


+ 1 








(y) 


0 










+1.6013 










+1.6013 


+ 1 






(z) 


0 


0 








- 1.2073 


.0 








+0.9578 


+0.5981 






-0.2495 


+0.5981 


+ 1 




1-90 1 


1.90 




0.5981 






.5263 




0.3148 




322.61 - 


192.96 = 


1 






322.61 0 . 


1 + 60.7438 = 


= 61.7438 






0. 1914 







Q 



'xz 



-0.2495 

-0.1313 



322.61 - 192.96 0^^ = 1.6013 

322.61 = 1.6013 - 25.3356 = - 23.7343 

^ = - 0.07357 

29.62 0^ - 47.43 + 35.76 0^^ = 1 

29.62 Qjjx " 1 ■ 3.4894 + 4.6953 = 2.2059 
0^ = 0.07447 



standard Errors of Unknovms after Adjustment: - 

= Ii n/ = 2.726 n/ 0.67447 = 2.726 (+ 0.2729) = + 0.7439 

\ ^ = 2.726 n/ 0.1914 = 2.726 (+ 0.4375) = + 1.1926 

\ = ds/ 0^2 = 2.726 n/ 0.5263 = 2.726 (+ 0.2294) = + 0.6253 




I 
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Results of least squares solution; - 
X = + 0.092899 = 

y = - 0.477200 = ^ 

z = - 1.131579 = 10,000 6f 
= + 0'.'09 + OV74 

5a = 1000 (- 0.4772) = - 477.2 m. 

' - O-O0O113158 

a^ = 6’, 378# 206 m. 

6a = -477 

New a = 6,377*729 m. + II90 m. 

f^ = 0.003390060 

8f = -0.0001131^8 

New f = 0.003276902 + 0.00006253 

l/f = 304.9 +5.4 

Tbe solution, as can be readily seen, does not conform to what 
is known to be the truth, and the standard errors resulting show that 
the precision of the solution is "somewhat" rovigh. This of course can 
be explained by noting that the length of the arc is quite short, and 
that the nximber of points used is small. Better results certainly 
could have been obtained if meiny points having greater longitude 
differences had been used, thus affording the benefit of the prime 
vertical components of the deflections. However, the method used is 



what is intended to be demonstrated here. 




t 
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It could have been expected that the results obtained from the 
data used would be as they are. Differences in latitude between the 
initial point £ind other points in the arc apparently differ systemat- 
ically. It is seen as a rotigh check of the standard error of the 
value of a obtained in the least squares solution, that if we take the 
average value of (cp^-cp) for the five northernmost stations and of the 
five southernmost stations, the resulting difference between the two 
is - 8'.'57« For the arc as short as the one used this is high. 

Taking the mean latitudes of each of the two groups of five 
points, we get = l6°26*. Taking the length of A (cp*-cp) at 31 

meters per second of arc, we have - 8'.'57 @ 31 m. = -266 meters. Using 
one radian as 57°3> we can determine the proportional difference for 
an arc of the length of the earth ^s radius ;- 

( - 266) = - 930 meters. 

This kind of rough calculation could tell us in advance 
approximately what results c£in be expected. 
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7.2 Dimensions from Sixteen Stations in South Central United States 

by the Astro-Gravimetric Method. 

For an illustration of the astro -gravimetric method as 
described in section 4.3# the published data on 16 points in South 
Central U. S. (including Kansas, Missouri, Oklahoma, Arkansas, Texas, 
and Louisiana) as computed by Rice (26), was used. 

In the referenced publication, all geodetic data were con- 
verted from Clarko'<^ 1866 Ellipsoid to the International Ellipsoid. In 
addition, as a result of the investigations. Rice determined that 
several of the astronomic observations which had been published in (7) 
were wrong, and were consequently done over. Moreover, all latitude 
and longitude determinations were reduced to the FK^ system of 
cataloging mean places of stars. The gravity anomalies used were free- 
air. 

All the data needed for solution of the present problem was 
given in Tables VI and VII of ref. (26), namely the latitudes and 
longitudes of each station, as well as the astro-geodetic and the 
gravimetric deflections for each. Here again, as in the problem of 
section 7»1, the extent of the area covered, and the number of stations 
involved are inadequate for a meaningful solution to the problem, but 
for illustrative purposes, the results are gratifying. 

lyfeade^s Ranch, Kansas, was again chosen as the initial point. 
The Vening Meinesz abbreviated formulas were determined to be adequate 
and were used for both | and q. Since Meade’s Ranch was one of the 16 
points given, the values of and are known, thus reducing the 



I 

I ' 



r 



I 



i 

i 
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mamber of unknowns to two (6N^ again being neglected). 

In adapting the formulas (^- 19 )> (^- 23 )> (^- 20 ), and (i(— 24 ) to 
this problem, the terms were grouped in the logical msinner for ease in 
computation of the coefficients as follow: - 

I = + d|^ + ( 9' " 9 ) - Ig 

I = + h^5l^ + (10,000) 6f + ( 9' - 9 ) - Ig 

= + cos (cp - cp^) 
h^ = + sin 9 sin (\-X^) 

= + [sin (9 - 9^) - 2 cos <p^ sin 9 sin^ ^ (4 - esc 1" 

[ 21 2 
[sin (9 - 9 q) - 2 cos cp^ sin 9 sin g- (\-\^)] sin 9^ 

- 4 cos 9 cos I (9 + 9 q) sin | (9 - 9^) 

- (2 . I tan <p^ Sin 4 <p„) sin (? - 

T) = dTljj + dTl^ + (\' - \) cos 9 g 

5 a 

= ggS ^l^’^o ^2 Too ^2 9 - ’Ig 

gg = - sin 9^ sin (X - X^) 

hg = + 1 

Jg = - cos 9^ sin (X- X^) ^ CSC 1 " 

2 

k_ = [ - cos cp Sin (X - X ) sin 9 
2 ^o ' o' ^o 

t J sin q.^ sin U <p^ sin X^) f] 
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Then the observation equations become: 



8 a 



5 = Ji ^ + ^1 (10,000) 8f + 1^ 



where = g, 8 | + h , 8 ^ + (| - | ) 

1 ^1 ^o 1 'o '^a 



= 9' - P 



55 =-(5 -5) =(5 -5) 

^o '^a ^g''o '^g ^a'^o 

&n = -(t1 -ti) =(q - l) 
o ' a 'g^o ' g a^o 



n = 02 ^2 + ^2 

where 1 ^ = 828 !^ + h^S 

T] = (X' - X) cos 9 

oL 




I 
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Table VI - Basic Data - Rice's l6 Points 



Point 

n 


Name 




1 


Meades R. 


39 


2 


Arcadia 


37 


3 


Dirks 


37 


4 


Twin 


36 


5 


Little Rock 


3k 


6 


Bogne 


33 


7 


Ecore 


33 


8 


Polk 


33 


9 


Bartley 


’ 33 


10 


Roby 


32 


11 


I>acasa 


32 


12 


Burns 


32 


13 


Sears 


32 


l4 


Brooks 


32 


15 


^num 


32 


16 


Legion 


32 



International Ellipsoid; - 



cp \ 



13' 


27" 


N. 98° 


32' 


31' 


38 


08 


94 


35 


47 


25 


4o 


99 


16 


21 


06 


57 


96 


46 


4o 


44 


51 


92 


16 


25 


k3 


05 


90 


53 


17 


36 


20 


92 


54 


25 


23 


05 


99 


05 


46 


12 


08 


101 


4o 


44 


44 


36 


100 


22 


46 


39 


05 


98 


4l 


30 


35 


28 


93 


14 


20 


33 


31 


100 


02 


16 


25 


09 


92 


07 


01 


19 


16 


101 


00 


00 


16 


19 


90 


06 


07 


a 


= 6, 


378,308 m. 







f = 1/29T 



= 0.003367003367 



Initial Point;- Meades Ranch 



cp^ = 39° 13' 26 '.'686 N. 
\ = 98° 32' 30'.’ 506 W. 



® = (5g - - ■ (5a - hk ‘ 

^ ’'0 * (''g - '“a'o = - - '>g)o = -(-°"30) = + 0L'30 



2 

k 

5 

6 

7 

8 

9 

10 

11 

12 

38 

14 

15 

16 

Poi] 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

l6 
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Table VI 


(cont'd) 










(1) 




(2) 


(3) 






(4) 


8l 




^1 




” ^g 








(0V90 g^) 


(0V30) 




(1) 


+ 


f f 

(2) + 


+ 


0.90 




0 


- 0.90 






C 


+ 


.90 


- 


0.01 


- .16 




+ 


0.'<3 


+ 


.90 


+ 


.00 


- .98 




- 


.C»d 


+ 


.90 


- 


.01 


- .73 




+ 


.16 


+ 


.90 


- 


.02 


+ .66 




+ 


1.54 


+ 


.90 


- 


.02 


+ .56 




+ 


1.44 


+ 


.90 


- 


,.02 


+ .07 




+ 


.95 




.90 


+ 


.00 


- .75 




+ 


.15 




.90 


+ 


.01 


- .89 




+ 


.02 


+ 


.89 


+ 


.01 


- .36 




+ 


.54 


+ 


.89 


+ 


.00 


- .25 




+ 


.64 


+ 


.89 


- 


.01 


+ .56 




+ 


1.44 


+ 


.89 


+ 


.00 


- .50 




+ 


.39 


+ 


.89 


- 


.02 


- .06 




+ 


.81 


+ 


.89 


+ 


.01 


- .45 




+ 


.45 


+ 


.89 


- 


.02 


+ .44 




+ 


1.31 


&2 


0 

<0 


^2 




T) - T) 

a 'g 






^2 


(0V90 gg) 


(o': 30 hg) 


If 


(1) 


t 

+ 


1 

(2) + 




0 


4* 


0.30 


0.30 






0 


+ 


0. oh 


+ 


.30 


- 1.24 




- 


0.90 


- 


.01 


+ 


.30 


- .34 




- 


.05 


+ 


.02 


+ 


.30 


- .61 




- 


.29 


+ 


.06 


+ 


.30 


0.00 




+ 


.36 


+ 


.08 


+ 


.30 


- .05 




+ 


.33 


+ 


.06 


+ 


.30 


+ '.76 




+ 


1.12 


- 


.01 


+ 


.30 


- .02 




+ 


.27 


- 


.03 


+ 


.30 


+ 2.26 




+ 


2.53 


- 


.02 


+ 


.30 


+ .74 




+ 


1.02 


- 


.00 


+ 


.30 


+ .85 




+ 


1.15 


+ 


.05 


+ 


.30 


+ .25 




+ 


.60 


- 


.01 


+ 


.30 


+ 1.14 




+ 


1.43 


+ 


.06 


+ 


.30 


+ .47 




+ 


.83 


- 


.02 


+ 


.30 


+ 1.88 




+ 


2.16 


+ 


.08 


+ 


.30 


+ .27 




+ 


.65 
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Table VII - Observation Eauations - Rice ' s l6 Points 



Point 








k 




1 


1 1 




0 




0 




0 


2 


• 


0.093 


+ 


0.672 


+ 


.73 


3 




.102 


+ 


.773 


- 


.08 


k 


- 


.176 


+ 


1.366 


+ 


.16 


5 


- 


,261 


+ 


1-993 


+ 


1.54 


6 


- 


.323 


+ 


2.484 


+ 


1.44 


7 


- 


.323 


+ 


2.554 


+ 


.95 


8 


- 


.329 


+ 


2.681 


+ 


.15 


9 


- 


.341 


+ 


2.767 


+ 


.02 


10 


- 


.366 


+ 


3.001 


+ 


.54 


11 


- 


.370 


+ 


3.050 


+ 


.64 


12 


- 


.379 


+ 


3.065 


+ 


1.44 


13 


- 


.376 


+ 


3.096 


+ 


.39 


14 


- 


.392 


+ 


3.146 


+ 


.81 


15 


- 


.390 


+ 


3.215 


+ 


.45 


16 


“ 


. 4 o 6 


+ 


3.205 


+ 


1.31 


n 1 




0 




0 




0 


2 


- 


0.172 


+ 


0 

• 

4 ^ 

-p- 

0 


- 


0.90 


3 


+ 


.032 


- 


.081 


- 


.05 


4 


- 


.077 


+ 


.197 


- 


.29 


5 


- 


.275 


+ 


.702 


+ 


.36 


6 




.334 


+ 


.850 


+ 


.33 


7 


• 


.246 


H 


.627 


+ 


1.12 


8 


+ 


.024 


- 


.062 


+ 


.27 


9 


+ 


.137 


- 


.349 


+ 


2.53 


10 


+ 


.080 


- 


.205 


+ 


1.02 


11 


+ 


.007 


- 


.017 


+ 


1.15 


12 


- 


.232 


+ 


.590 


+ 


.60 


13 




.065 


- 


,167 


+ 


1.43 


14 


- 


.28c 


+ 


.714 


+ 


.83 


15 


+ 


.107 


- 


.274 


+ 


2.16 


16 


- 


.368 


+ 


.937 


+ 


.65 




f 
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Taole VIII 


lonot Squaroc 


Solution 


• Sice ' s 


16 Points 


6a 

100 


10,000 6f 


1 


L 


1 (a) 1(f) 


6a 

100 + 2.17 
1 


- 14.21 
+ 0.54839 


- 3.62 
+1.66820 


- 15.66 
+7.21659 


1 


10,000 6f 


+ 106.07 
- 93.05 


+ 27.95 
- 23.71 


+ 119.81 
- 102.55 


0 

+6.548 




+ 13.02 
1 


+ 4 . 2 '!- 

-0.32565 


+ 17.26 
-1.32565 


+6.548 1 


1 




+ 30.84 

- 6 .o 4 

- 1.38 


+ 55.17 

- 26.12 

- 5.62 




[w] 




+ 23.42 


+ 23.43 







10,000 6 f 


Check 


6 a 

100 


Check 


1 


-0.32565 


-1.32565 


+1.66820 


+7.21639 


10,000 6 f 






- 2.13248 


-8.68087 


Z 


-0.32565 


-1.32565 


- 0.46428 


- 1.46426 



Standard Error of Single Observations- 

- °- 5°3 

(+ 2.17 Q - (0.503) - 1 

„ 8.148 

' 2.755 



= 0 



Srf 13.02 " 0.077 

standard Errors of Unknowns j- 



m 



= u n/ Q 



m 

= 0.8836 = 1.712, 

m = 171.2 m 

3 « 



10,000 m^=0. 8836 n/o. 077=0. 24520 
= 0.000024520 



68 



Results of Least Squares Solution 



10,000 5 f = - 0.32565 


+ 0.24520 


6 f = - 0.000032565 


+ 0.000024520 


f = 0.003367003 

0 




f = f + 6 f = 0.00333^438 
0 


+ 0.000024520 


J = 299.9 + 2.2 




6 a 




100 = - 0.46428 + 


1.71 meters 



5 a = - 1<6 + 171 

= 6,378,388 

a = a^ + 6 a = 6,378,3^2 + I7I meters 

Here, as in the problem of section 7*1, the standard errors 
of the corrections are almost as large as or larger than the 
corrections themselves. But the fact that they, as well as the 
corrections to f and a themselves, are relatively small, testifies as 
to the value of the method. However, it must be realized that the 
standard errors reflect only the precision of the least squares 
adjustment, asstiming the given data has no error, which of course is 



not true. 
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7.3 Computation of the Flattening from Gravity Anomalies. 

Given a set of isostatic gravity anomalies (Hayford, 

D = 113 «T km.) at 67 points at various random positions scattered 
over the earth, and computed on the basis of Helmert's gravity 
formula of I90I, we can demonstrate the solution of a problem of 
correcting the gravity formula and of obtaining a value for the 
flattening of the ellipsoid (13 )> p. (l^)* The method is that 

of Chapter tut extended to ..include the longitude term of the 
formula. 

Helmert's Formula? - 

7 = 978.030 (1 + 0.005302 sin^cp + 0.000007 sin^2cp) 

for which the p term yields the value of f = I/298.2. See 

(13), p. 78. 

The general form of the gravity formula, including the 
longitude term is;- 

2 2 2 

7 = 7„[l + P sin cp + € sin 2cp + r cos cp cos 2(\ - \ )] 
jlj o 

2 2 2 

= 7g[l + P sin cp+esin 2qp+r cos cp (cos2\cos2\^+ sin 2\ sin 2\^)] 

Making corrections, we have 7' = 7 + A7 , 

2 2 2 

where = x' + y'sin cp + z'cos cp cos 2\ + u'cos qp sin 2\ 
x' = tmknown corr'n to 7^, in mgals 
y' = " " " P in mgals («■ 1.022 x lO"^) 

z' = " " " r cos 2\^ in mgals (» 1.022 x lo"^) 

u' = " " " r sin 2\^ in mgals 1.022 x lO’^) 



TO 



Using = - V, we have 

2 2 2 
x' + y'sin cp + z'cos cp cos 2 \ + u’cos cp sin 2 \ 

The observation equations will then bes- 

ax* + by' + CZ* + du' - l= -v 

where a = 1 



^ . 2 
b = sin <p 

2 

c = cos cp cos 2 \ 
2 

d = cos cp sin 2 X 

1 = ^ 



Zg = - V 



The solution, according to the method of least squares will 

2 

be £ V = minimum. 



* 



% 



* i I • 






Poi: 

1 

2 

3 

k 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

l6 

17 

18 

19 

20 

21 

22 

23 

2 k 

25 

26 

27 

28 

29 

30 

31 

32 

33 
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- Isostatlc Gravity Anomalies (Rayford, D = 113« 



<p X 



79?5 


11?5 


+ 16 


79.5 


17.5 


+ 8 


76.5 


15.5 


- 26 


71.5 


27.5 


+ 20 


65.5 


?4.5 


- 5 


64.5 


11.5 


+ 28 


62.5 


17.5 


- 20 


60.5 


22.5 


+ 24 


58.5 


38.5 


+» 16 


57.5 


40.5 


+ 28 


55.5 


13.5 


+ 10 


54.5 


39-5 


- 11 


54.5 


55.5 


+ 4 


53.5 


29.5 


+ 15 


53.5 


50.5 


+ 23 


52.5 


25.5 


+ 15 


>1-5 


- 0.5 


+ 9 


^1.5 


11.5 


+ 36 


50.5 


10.5 


+ 32 


48.5 


16.5 


+ 31 


47.5 


39-5 


+ 43 


46.5 


7.5 


+ 20 


46.5 


34.5 


+ 6 


45.5 


10.5 


+ 30 


44.5 


- 1.5 


+ 6 


44.5 


39.5 


+ 31 


43.5 


42.5 


+ 43 


40.5 


44.5 


+ 74 


39.5 


53.5 


+ 15 


36.5 


+ 6.5 


+ 14 


35.5 


14.5 


+ 70 


28.5 


33.5 


+ 17 


26.5 


33.5 


+ 31 



Point 


(p 


X 


34 ' 


22?5 


36?5 


35 


18.5 


38.5 


36 


14.5 


- 17.5 


37 


13.5 


42.5 


38 


- 15.5 


38.5 


39 


- 26.5 


15.5 


4o 


50.5 


120.5 


4l 


49.5 


94.5 


42 


48.5 


116.5 


k 3 


47.5 


96.5 


44 


47.5 


79.5 


45 


46.5 


123.5 


46 


46.5 


92.5 


47 


46.5 


60.5 


48 


45.5 


102.5 


49 


44.5 


110.5 


50 


44.5 


74.5 


51 


43.5 


94.5 


52 


43.5 


74.5 


53 


42.5 


71.5 


54 


39.5 


120.5 


55 


39.5 


84.5 


56 


38.5 


78.5 


57 


36.5 


87.5 


58 


35.5 


105.5 


59 


34.5 


79.5 


60 


31.5 


94.5 


61 


56.5 


59.5 


62 


51.5 


104.5 


63 


46.5 


61.5 


64 


38.5 


70.5 


65 


23.5 


79.5 


66 


22.5 


77.5 


67 


11.5 


78.5 



km. ) 






+ 26 
+ k & 
+ 44 

+ 43 
+ 27 
+ 34 
- 1 

+ 9 

+ 10 

+ 19 

+ 6 

- 5 

+ 58 

- 1 
+ 23 

+ 14 
+ 26 
+ 14 

- 16 
+ 1 

- 20 
- 11 

- 5 
+ 14 
+ 11 
+ 10 

- 4 
+ 55 

- 2 
+ 14 
- 66 
+ 18 
+ 27 
- 4o 
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Table X - Error Equations - Flattening and Correction 
to Gravity Formula frgm Gravity Anomalies 



3 int 


a 


b 


c 


d 


1 




(X-) 


(y) 


(z’) 


(u' ) 


= - ^ 


1 


1 


.9669 


.0306 


.0130 


- 16 


2 


1 


.9669 


-0272 


.0190 


- 8 


3 


1 


.9^56 


.0467 


.0281 


+ 26 


k 


1 


.8993 


.o «^77 


.0825 


- 20 


5 


1 


. 8281 ’ 


.1128 


.1298 


+ 5 


6 


1 


.8147 


.1706 


.0724 


- 28 


7 


1 


.7868 


.1746 


.1223 


+ 20 


8 


1 


.7576 


.1714 


.1714 


- 24 


9 


1 


.7269 


.0614 


.2660 


- 16 


10 


1 


.7113 


.0452 


.2851 


- 28 


11 


1 


.6791 


.2858 


.1456 


- 10 


12 


1 


.6628 


.0643 


.3310 


+ 11 


13 


1 


.6628 


- .1209 


.3148 


- 4 


Ik 


1 


.6463 


.1822 


.3033 


- 15 


15 


1 


.6463 


- .0062 


.3537 


- 23 


l 6 


1 


.6295 


.2332 


.2880 


- 15 


17 


1 


.6125 


.3874 


- . 00 ^^ 


- 9 


18 


1 


.6125 


.3567 


.1514 


- 36 


19 


1 


.5954 


.3778 


.1450 


- 32 


20 


1 


.5610 


.3682 


.2391 


- 31 


21 


1 


.5436 


.0871 


.4480 


- 43 


22 


1 


.5262 


.4577 


.1226 


- 20 


23 


1 


.5262 


.1698 


.4424 


- 6 


2 k 


1 


.5088 


.4586 


.1761 


- 30 


25 


1 


.4913 


.5081 


- .0266 


- 6 


26 


1 


.4913 


.0971 


.4994 


- 31 


27 


1 


.4739 


.0459 


.5242 


- 43 


28 


1 


.4217 


.0101 


.5781 


- 74 


29 


1 


. 4 o 46 


- .0726 


.5909 


- 15 


30 


1 


.3538 


.6297 


.1454 


- 14 


31 


1 


.3372 


.5796 


.3213 


- 70 


32 


1 


.2277 


.3017 


.7109 


- 17 


33 


1 


.1991 


.3129 


.7372 


- 31 




-nsjf. 








>int 

34 

35 

36 

37 

38 

39 

4o 

4i 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

bO 

61 

62 

63 

64 

65 

66 

67 
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a 


b 


(x*) 


(y) 


1 


.1465 


1 


.1007 


1 


.0627 


1 


.0545 


1 


. 0714 . 


1 


.1991 


1 


.5954 


1 


.5782 


1 


.5610 


1 


.5436 


1 


.5436 


1 


.5262 


1 


.5262 


1 


.5262 


1 


.5088 


1 


.4913 


1 


.4913 


1 


.4739 


1 


.4739 


1 


.4564 


1 


.4046 


1 


. 4 o 46 


1 


.3875 


1 


.3538 


1 


.3372 


1 


.3208 


1 


.2730 


1 


.6954 


1 


.6125 


1 


.5262 


1 


.3875 


1 


.1590 


1 


.1465 


1 


.0398 



X (cont*d) 



c 

(z' ) 



.2496 

.2023 

.7678 

.0825 

.2089 

.6865 

- . 1962 

- . 4165 

- .2642 

- .4448 

- .4261 

- .1852 

- . 4721 

- .2441 

- .4452 

- .3840 

- .4361 

- .5197 

- .4511 

- .4341 

- .2886 

- .5844 

- .5638 

- .6438 

- .5681 

- .6340 

- .7I8C 

- .147, 

- .338s 

- .2581 

- .4759 

- .7852 

- .7736 

- .8839 



d 1 



(u*) 




.8163 


- 26 


.8763 


- 48 


.5376 


- 44 


.9420 


- 43 


.9048 


- 27 


.4124 


- 34 


.3539 


+ 1 


.0660 


- 9 


.3506 


- 10 


.1027 


- 19 


.1636 


- 6 


.4362 


+ 5 


.0413 


- 58 


.4062 


+ 1 


.2076 


- 23 


.3338 


- 11 


.2620 


- 26 


.0823 


- 14 


.2710 


+ 16 


.3271 


- 1 


.5207 


+ 20 


.1136 


+ 11 


.2393 


+ 5 


.0564 


- 14 


.3413 


- 11 


.2434 


- 10 


.1137 


+ 4 


.2664 


- 55 


.1879 


+ 2 


.3975 


- 14 


.3854 


+ 66 


.3014 


- 18 


.3607 


- 27 


.3752 


+ 4o 
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Table XI 



Least Squares Solution 



Flattening and Correction 



to Gravity Formula from Gravity Anomalies 



Forward Solutions- 





x‘ 


y’ 


z* 


u* 




1 




z: 


X* 


+67 


+ 33.20 


- 4.17 


+ 13.17 




1033 




923.80 




- 1 


-0.4955 


+0.0622 


-0.1966 


+ 


15.4179 


+ 


13.7881 


y' 




+ 19.76 


- 1.29 


+ 4.87 




460.99 




4o4.45 






- 16.45 


+ 2.07 


- 6.53 


+ 


511.85 


+ 


457.74 


(1) 




+ 3.31 


+ 0.78 


- 1.66 


+ 


50.86 


+ 


53.29 






- 1 


-0.2356 


+0. 5015 


- 


15.3656 


- 


16.0997 


z* 






+ 10.73 


+ 0.94 




183.53 




177.32 








- 0.26 


+ 0.82 


- 


64.25 


- 


57.46 








- 0.18 


+ 0.39 


- 


11.98 


- 


12.55 


(2) 






+ 10.29 


+ 2.15 




259.76 




247.33 








1 


-0.2089 


+ 


25.2439 


+ 


24.0350 


u* 








+ 9.65 


• 


337.37 


• 


308.74 










- 2.59 


+ 203.09 


+ 


181.62 










- 0.83 


+ 


25.51 


+ 


26.72 










- 0.45 


+ 


54.26 


+ 


51.67 


(3) 








+ 5.78 




54.51 




48.73 










1 


+ 


9.4308 


+ 


8.4308 



+52,563 

-15,926.69 

- 781. 49 

- 6557.36 

- 514. 07 



+50,548.11 

-14,243.06 

- 818.83 

- 6243.32 

459.56 



[w] 



+28,783.39 +28,783.34 



Rack Solutions- 
' Check 



u 



Check 



y’ 



Check 



Check 



1 +9.4308 +8.4308 +25.2439 +24.0350 -15.3656 -16.0997 +15.4179 +13.7881 

u* - 1.9701 - 1.7612 + 4.7295 + 4.2280 - 1.8541 - 1.6575 

z» - 5.4833 - 5.2477 ♦ 1.4k76 + 1.3854 

’ + r.9872 + 8.4827 

+9.4308 +8.4308 +23.2738 +22.2738 -16.1194 -17.1194 +22.9986 +21.9987 



Table XII - Correlate Equations - Flattening> etc.> f rom 
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• rH 
OO 
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s 

O O O CVJ 



0 
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iT\ iTv 

oo on 

O O CVJ CVJ 



0 

1 



o 

9 



O -=^ 
O O iPv O 

9 • 

o o 
+ + 



X 



ir\ 

O 

0 

1 



iC\ 

0 

1 



CVJ 

CVJVQ 

rH 

> O rH 

• • 

o o 
+ + 



o\ 

00 

rH 

o 

+ 



VO if\^ 
VO cb 
oS-^ oo 

O H CVJ O 

• • e 

0 o o 

1 I I 



N 



CVJ 



oo 



X 



5.78 - 0.4825 +O.55OT -0.2089 



Solution of Correlate Numbers 



5.78 = - 0.4825 

0^ = - 0.08478 

10.29 Q^z " = 0.1789 



= (0.1789 + 0.1823) ♦ 10.29 = 0.03510 
3.31 + 0.78 (.03510) - 1.66 (- .08478) = - 0.4955 

0^ = (- 0.4955 - 0.0274 - 0 . 14 o 7 ) * 3.31 = - 0.20048 
67 Q + 33.20 (- .20048) - 4.17 (.03510) + 13.17 (- .08478) 
0^ = (1 + 6.6559 + 0.1464 + 1.1166) ♦ 67 = 0.13312 
5.78 = 0.5507 

Qy^ = 0.09528 

10.29 0^2 (.09528) = - 0.2356 

Qy^ = (- 0.2356 - 0.2049) ♦ 10.29 = - 0.04281 
3.31 0^ + 0.78 (- .04281) - 1.66 (.09528) = 1 
Qyy. = (1 + 0.0334 + 0.1582) ♦ 3.31 = 0.36000 



5.78 Q, 



zu 



0.2089 



Q = - 0.03614 

zu 

10.29 + 2.15 (- .03614) = 1 



0^^ = (1 + 0.0777) ♦ 10.29 = 0.10473 

5.Y8 V = 1 

o - 0.17301 



;• « 




I 



-p. 






" i 










standard Error of an observation of unit weight: - 



^ = vT 4507 = + 21.37 mgal 



Standard Errors of the unknowns after adjustment; - 
= PnT , etc. 

= 21.37 -J 0.1331 = + 7.78 (mgal. ) 

m ' = 21.37 n/ 0.3600 = + 12.82 (mgal. * 1.022 x 10 

y 

m ' = 21.37 v/ 0.1047 = + 6.92 (mgal. ♦ 1.022 x lO"^) 

m^' = 21.37 n/ 0.1730 = '+ 8.89 (mgal. » 1.022 x lo'^) 

The values for the unknowns, solved in Table XI, and their 
standard errors, are here converted to gals.:- 



X 


= X* 


X lo"^ 


y 


= y' 


X (1.022) X 10"^ 


z 


= z' 


X (1.022) X lO"^ 


u 


= u* 


X (1.022) X 10"^ 



= + 0.0230 + 0.0078 gals. 
= - 0.000016+0.000013 gals. 
= + 0.0000238+ 0.0000071 gals. 
= + 0.0000096+ 0.0000091 gals. 



The components of the longitude term are:- 

r^ = + u^ = 0.0006586 X 10"^ > r = 0.0000257 

2 X = 22°, X = 11° 
o ^ o 



Applying the computed corrections : - 

= 978.03c + X = 978.053 gals. 

S = 0.005302 + y = 0.005286 
€ 



= 0.000007 



(unchanged) 
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Sind addixig the longitude term;- 

0.000026 cos^cp cos 2 (\ - 11°) 

we obtain the corrected gravity formula: - 

r = 978.053 [1 + 0.005286 sin^cp + 0.000007 sin^ 2 9 
+ 0.000026 cos^9 cos 2 (\ - 11°)]. 



Next, the flattening f is solved from the corrected value of p. 

At this point it should be pointed out that if the principal problem 
is to find the flattening, it’would not be necessary to include 
solution of the longitude term. With slight difference the same 
result for the correction to p is obtained when only two unknowns exist 
in the normal equations. Witness, using the first reduced normsQ. 
equation in Table XI:- 

3*31 y’ + 50.86 = 0 

y* = - 50.86 * 3.31 = - 

X = - 15 A X 1.022 X lO"^ = - 0.000016 

P = 0.005302 - 0.000016 = 0.005286 



Now using Clairaut's Formula, (5-1 )> we obtain f fron the new 
value of p. 

0.0052855 = |m-f-^mf 



m 



m 



2 

0) a 



2 7 T -1 

= mm 



_ f 2 IT .2 637^820,000 cm. 

■ ^86164'' 978.05 galT 

39.^78 X 637.82 _ 25,179.858 

"" 7424.2 X 978.05 7 , 261,124 " 



0. 003467 b 



I 









] 
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0. 0052855 

f 



l/f 

0 

i"f 

Of course, the result reflects the fact that only a few values 
of ^ were used, and at scattered locations. The solution by 
Heiskanen, referred to in Chapter 5> shows that the use of a great many 
points results in a value of f equal to that obtained from Helmert’s 
original formula. The results in the computation given here are sur- 
prising considering the fact that the points chosen were completely 
random with no thought given to even approximately even distribution. 

7.4 Computation of Gravimetric Deflections 
7.41 Procedure sind General Information. 

The procedure used for computation of the gravimetric deflect- 
ions at the two points chosen, is that involving the use of the Rice 
Circle -Rings. The extent of readily available gravity material 
dictated that the area covered in both cases was to be withir. the 
circle-ring no. 49, having a radius of 399.0 km. 



I (0.0034678) - f [1 + ^ (0.0034678)] 
0.008670 - f (1 + 0.0042109) 
0.008670 - 0.005286 



1.0042109 

0.0033698 + 
= 296.75 + 0.74 



0.0033840 

1.0042109 



0.0000083 



0.0052855 + 0.0000130 
= + “1^ (0.0000131) = + 0.0000083) 
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A template was prepared to the scale of the available anomaly 
maps, having circles of radii listed in Table I and with 10° sectors. 

The usual manner of orienting the template, which is to have zero 
azimuth bisect one of the sectors making mean azimuths of all sectors 
multiples of 10°, was not followed here because the writer’ s method of 
assuring proper sign determination dictated that each compartment have 
its own quadrant designation. Therefore 1±ie azimuth lines were 
nvimbered in multiples of 10°, and the mean sector azimuths are multi- 
ples of 

For the point in the northern hemisphere, zero azimuth is nortxj, 
and for the one in the southern hemisphere it is south. Following the 
usual practice for points in Europe, \ and t) were considered positive 
eastward, so for both points, the azimuths were measured eastward from 
zero| thus they are measured clockwise in the northern hemisphere, and 
counterclockwise in the southern hemisphere. This will preserve signs 
automatically, using the same procedures for both. Figure 12 illu- 
strates the azimuth convention used. 

Remembering that the plumb line (below the surface) is pulled 
toward the mass surplus and thus in the direction of the more positive 
of the anomalies on opposite sides of the computation point, we can 
devise a system of assuring correct sign. Figure 13 illustrates this. 

I is considered positive toward the nearest pole, tj is considered 
positive eastward. Table XIII gives criteria to be followed. 
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//0/ZTHEJ^A/ 

HEMfSPHE/EE 



Southern 

HEM/SPHBR£. 



Fig. 12 



Azimuth Convention for Computing 
Gravimetric Deflections for the 
Points Computed in this Thesis. 
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A/e£iresi~ PoLe. 




, 2ir ^ TT I ^ 




(A^j. -A^^) A. O 



Fig. 13 
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Table XIII - Criteria for Determining Effective Azimuth 

and Sign of Deflection 



Quadrants 


Anomaly Jradient 


Effective 

Azimuth 


cosine 

and 

1 


sine 

and 


I and III 


(agl -:5gjjj)<0 




> 0 


> 0 




(igj - Agjjj) > 0 


^:II 


< 0 


< 0 


II and IV 


(^JJ - Agj^) < 0 




< 0 


> 0 




(^II - > 0 




> 0 


< 0 

„ -- 



Computation proceeds, beginning in the first quadrant with the 
effect of that quadrant and the one opposite, the third. After com- 
pleting the first and third, the second quadrant is begun, determining 
the effect of that quadrant and the one opposite it, the fourth. Thus 
the entire 3^0° in azimuth are covered. 

In order to obtain data for a graph of the effects of areas of 
different radii on the deflection of a point, the computation was made 
grouping the effects of four zones together from the point outward. 

Then accumulative effects were computed and plotted for both | and q. 

The inner circle, for one case 5*125 km. in radius, and for the 
other 10.15 km., was computed by the one-gradient method. Insufficient 
local gravity data at the points dictated that the three -gradient 
method would yield results of meaningless accuracy. 

The two points conputed were Laiska, Finland (64° 03' N, 

28° 50' E), and Wonderfontein, Union of South Africa (25° 46' S, 
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29° 53' E). The points were chosen to contribute data for the com- 
putation of the equatorial radius from the arc of the 30th meridian. 

In addition, their locations had to be such that the anomaly maps 
available to the writer would provide data to a large enough radius 
aroxmd each point to make computation by the circle -ring method 
feasible. The result was that both the points turned out to be scme- 
what removed from areas of concentrated gravity stations. However, 
point-anomaly values obtained from files of the Institute of Geodesy, 
Photogrammetry, and Cartography, The Ohio State University, for the 
general area around the points, allowed the plotting of iso -anomaly 
contour maps to the scale 1" = 6 km. (l ; 240, OOO) for both points, and 
the writer proceeded with the computation \mdaunted by the lack of 
enough gravity data. 

The small scale maps used for the bulk of the computation were 
the following: - 

Lai ska : Free -air anomaly map of Finland on the scale 

1 : 1,000,000, covering the area 21° - 35 ° East longitude 
(extending into Russia) and 6o° - 69° North latitude. 

Contours for the Gulf of Bothnia were sketched in from a 
smaller scale map of that area. The maps used were 
prepared from gravity determinations of the Finnish 
Geodetic Institute and were computed on the International 
Gravity Formula. Contour interval was 5 mgal. 

Wonderfontein ; Bouguer anomaly lithographed map sheets of 
Union of South Africa, Northeast and Southeast sheets, on 
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scale 1 : 1,000,000, 10 mgal contours; total coverage of 
the two sheets i 22° 00' to 30' S, 23° 30' to 33° 50' 

E. These maps were published in South Africa. 

Topographic maps published as World Aeronautical 
Charts by U. S. Air Force (ACIC), scale 1 ; 1,000,000; 
contour interval is 1000 ft. , but adequate point elevations 
to the nearest foot are given for the use to which the maps 
were put. Total coverage of sheet numbers 1275, 1276, 1299# 
1300, 1397, and 1398 is about the same as that of the 
Bouguer anomaly maps. 



Since the South African point was estimated from a Bouguer 
anomaly map, it was necessary to use the topographic maps also, for the 
purpose of obtaining mean elevations for each compartment of the same 
Rice ring template used for the anomalies. Thus armed with the 
elevations, the Bouguer correction could be computed and subtracted 
from the Bouguer anomaly to obtain the Free Air anomaly for the com- 
partment. 



= Zigg - g'x, = Free Air anomaly 



B 



where 

g' 



B 



B 



Bouguer anomaly 
Bouguer correction 

I'" 



m 



= density of cr^ist assxamed for area 

= 2.67 



P 



m 



mean density of earth 

5.576 
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S'b = | - ; 5 y6 (o* 09^06 ) h mgals 

= - 0.03378 h mgals. 



(h in feet) 



7.^2 Computation for Laiska, Finland. 

The computations for the circle-ring zones, proceeding from the 
outside zone, no. 48, and continuing inwaixi toward the computation point 
to zone no. 27, are included in Table XIV. 

The inner circle, r^ = 10.15 toi. , is computed below. The 
effects of both the inner circle and the circle-ring zones are 
accumulated in Table XV. Curves of the effect on deflections of in — 
creasing gravity field radius are shown in figure 15 . 




Fig. l4 - Inner Circle - Laiska. 



Using the values shown in fig. l4, we obtain: - 



d|" = 0.0525 (^3 - 



= 0.0525 (2.2 - 2.0) = 0.0525 (0.2) = + O'.'OlO 

dT)'' = 0.0525 (^ - -^gg) 



(- 0.3 - 4.8) = 0.0525 (- 5.1) = - OV268 



0.0525 
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Table XIV - Computation of Gravimetric Deflections 

at Lai ska ^ Finland 



^ (mgal) 



Quadrant 




T 


III 


Sector 

Zone 


0-10 


180 -19c 


48 




- 10 


- 30 


47 




- 20 


- 5 


E Zig 




- 30 


- 35 


I - III 






+ 5 


d© @ 0"001 






0'.' 005 


Eff. A 






185 


cos A 






- .9962 


Bin A 






- .0872 


H 






- 0"005 


dT) 






o'.'ooo 


46 




- 15 


+ 15 


45 




0 


+ 25 


44 




+ 10 


+ 30 


43 




+ 10 


+ 20 


E 




+ 5 


+ 90 


I - III 




- 85 




d© 




0.085 




A 




5 




cos 




.9962 




sin 


+ 


.0872 




d| 


+ 


.085 




dn 


+ 


.007 




42 




+ 15 


+ 15 


4 l 




+ 18 


+ 11 


4 o 




+ 18 


+ 10 


39 




+ 15 


+ 8 


E 




+ 66 


+ 44 


I - III 






+ 22 


d© 






0.022 


A 






185 


cos 






- .9962 


sin 






- .0872 


H 






- .022 


dn 






- .002 



I 


III 


I 


10 - 20 


190-200 


20 - 30 


- 15 


- 30 


- 20 


- 25 


- 15 


- 10 


- 4 o 


- 45 


- 30 




+ 5 


- 10 




0"005 


OVOlO 




195 


25 




- .9659 


+ .9063 




- .2588 


+ .4226 




- 0"005 


+ o'.’oo9 




- o'iool 


+ 0'.'004 


- 10 


+ 20 


0 


+ 5 


+ 20 


+ 5 


+ 10 


+ 15 


+ 10 


+ 10 


+ 15 


+ 15 


+ 15 


+ 70 


+ 30 


- 55 




- 30 


0.055 




0.030 


15 




25 


+ .9659 




+ .9063 


+ .2588 




+ .4226 


+ .053 




+ .027 


+ .014 




+ .013 


+ 15 


+ 15 


+ 18 


+ 20 


+ 13 


+ 22 


+ 20 


+ 10 


+ 25 


+ 20 


+ 5 


+ 24 


+ 75 


+ 43 


+ 89 




+ 32 






0.032 






195 






- .9659 






- .2588 






- .031 






- .008 





III 

200-210 



- 20 

0 

- 20 



+ 20 
+ 20 
+ 10 
+ 10 

+ 6 o 



+ 7 

+ 10 
+ li^ 

+ 7 

+ 38 

+ 51 
0.051 
205 

- .9063 

- 1+226 

- ,.oli 6 

- .022 




r * 
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E 

I - 

de 

A 

cos 

sin 

d| 

dTi 



E 

I - 

de 

A 

cos 

sin 

dTi 



E 

I - 

de 

A 

cos 

sin 

dTi 









Table 


XIV 


(cont'd) 














I 


III 


I 




III 


I 






Ill 




30 - 4 o 


210-220 


ko - 


50 


220-230 


50 - 


60 


230-240 


48 




+ 5 


+ 10 


+ 


20 


0 


+ 


35 




- 10 


47 




+ 10 


+ 10 


+ 


30 


- 5 


+ 


45 




- 10 






+ 15 


+ 20 


+ 


50 


- 5 


+ 


80 




- 20 


III 




- 5 








+ 55 








+ 100 






.005 








.055 








.100 






35 








225 








235 




+ 


.8192 








- .7071 






- 


.5736 




+ 


.5736 








- .7071 






- 


.8192 




+ 


.004 








- .039 






- 


.057 




+ 


.003 








- .039 








.082 


46 




+ 5 


+ 5 


+ 


20 


+ 5 


+ 


25 




+ 10 


45 




+ 5 


+ 5 


+ 


10 


+ 10 


+ 


15 




+ 5 


44 




+ 10 


+ 10 


+ 


10 


0 


+ 


10 




0 


43 




+ 15 


+ 5 


+ 


10 


- 5 


+ 


10 




0 






+ 35 


+ 25 


+ 


50 


+ 10 


+ 


60 




+ 15 


III 






+ 10 






+ 4 o 








+ 45 








.010 






.o 4 o 








.045 








215 






225 








235 








- .8192 






- .7071 






- 


.5736 








- .5736 






- .7071 






- 


.8192 








- .008 






- .028 








.026 








- .006 






- .028 






- 


.037 


42 




+ 18 


+ 6 


4 - 


17 


+ 5 


+ 


14 




+ 5 


41 




+ 22 


+ 12 


+ 


22 


+ 12 


+ 


17 




+ 10 


4 o 




+ 27 


+ 15 


+ 


26 


+ 15 


+ 


20 




+ 5 


39 




+ 34 


+ 13 


+ 


28 


+ 13 


+ 


20 




+ 6 






+ 101 


+ 46 


+ 


93 


+ 45 


+ 


71 




+ 26 


III 






+ 55 






+ 48 








+ 45 








.055 






.048 








.045 








215 






225 








235 








- .8192 






- .7071 






- 


.5736 








- .5736 






- .7071 






- 


.8192 








- .045 






- .034 






- 


.026 








- .032 






- .034 






- 


.037 




► Itn 
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Table XIV (cont'd) 



I III 







6o - 


70 


240-250 




48 


+ 


45 


- 15 




47 


+ 


40 


- 15 


z 




+ 


85 


- 30 


I - 


III 






+ 115 


de 








.115 


A 








245 


cos 








- .4226 


sin 








- .9063 










- .049 


dTl 








- .104 




48 


+ 


15 


- 10 




45 


+ 


10 


0 




44 


+ 


10 


+ 5 




43 


+ 


10 


- 5 


Z 




+ 


45 


- 10 


I - 


III 






+ 55 


d0 








.055 


A 








245 


cos 








- .4226 


sin 








- .9063 


U 








- .023 


dTl 








- .050 




42 


+ 


13 


+ 5 




4l 


+ 


12 


+ 5 




4o 


+ 


10 


0 




39 


+ 


10 


0 



I 




III 


I 




III 


D - 


80 


250-260 


80 - 


90 


260-270 


+ 


45 


- 4o 


+ 


35 


- 4c 


+ 


35 


- 25 


+ 


30 


- 35 


+ 


80 


- 65 


+ 


65 


- 75 






+ 145 






+ l4o 






.145 






.140 






255 






265 






- .2588 






- .0872 






- .9659 






- .9962 






- .038 






- .012 






- .l4o 






- .139 


+ 


15 


- 20 


+ 


20 


- 25 


+ 


10 


- 15 


+ 


15 


- 25 


+ 


10 


- 10 


+ 


20 


- 10 


+ 


15 


- 5 


+ 


20 


- 5 


+ 


50 


- 50 


+ 


75 


- 65 






+ 100 






+ l4o 






.100 






.140 






255 






265 






- .2588 






- .0872 






- .9659 






- .9962 






- .026 






- .012 






- .097 






- .139 


+ 


17 


0 


+ 


15 


- 3 


+ 


10 


+ 3 


+ 


5 


0 


+ 


5 


0 


+ 


4 


0 




6 


0 


+ 


7 


0 



Z 


+ 45 +10 


+ 38 +3 


+ 31 - 3 


I - III 


+ 35 


+ 35 


+ 34 


d 0 


.035 


.035 


.034 


A 


245 


255 


265 


cos 


- .4226 


- .2588 


- .0872 


sin 


- .9063 


- .9659 


- .9962 


d| 


- .015 


- .009 


- .003 


dTl 


- .032 


- .034 


- .034 
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Table XIV ( c ont ' d ) 



II IV 





90 - 


100 


270-280 


k & 


+ 


30 


- 15 


hi 


+ 


20 


- 30 


z 


+ 


50 


- 45 


II - IV 






+ 95 


de 






.095 


A 






275 


cos 






+ .0872 


sin 






- .9962 


d| 






+ .008 


dll 






- .095 


h 6 


+ 


15 


- 30 


45 


+ 


10 


- 25 


44 


+ 


15 


- 5 


43 


+ 


10 


+ 5 


E 




50 


- 55 


II - IV 






+ 105 


d 0 






.105 


A 






275 


cos 






+ .0872 


sin 






- .9962 


d| 






+ .009 


dT) 






- .105 


42 




0 


2 


4 l 


+ 


7 


- 2 


4 o 


+ 


15 


- 2 


39 


+ 


15 


1 


E 


+ 


37 


- 7 


II - IV 






+ 44 


d 0 






.044 


A 






275 


cos 






+ .0872 


sin 






- .9962 


d| 






+ .oo 4 


dTi 






- .044 



II 




IV 


II 




IV 


100- 


110 


280-290 


110- 


120 


290-300 


+ 


15 


- 20 


+ 


10 


- 20 


+ 


10 


- 30 


+ 


5 


- 35 


+ 


25 


- 50 


+ 


15 


- 1^)5 






+ 75 






+ 70 






.075 






.070 






285 






295 






+ .2588 






+ .4226 






- .9659 






- .9063 






+ .019 






+ .030 






- .072 






- .063 


+ 


5 


- 35 




0 


- 25 




0 


- 25 


- 


5 


- 25 




0 


- 5 


- 


10 


- 20 


- 


5 


+ 5 


- 


5 


- 15 




0 


- 60 




20 


- 85 






+ 60 






+ 65 






.060 






.065 






285 






295 






+ .2588 






+ .4226 






- .9659 






- .9063 






+ .016 






+ .027 






- .058 






- .059 




0 


- 7 




0 


- 5 


+ 


12 


- 7 


+ 


12 


- 8 


+ 


20 


- 6 


+ 


18 


- 10 


+ 


20 


- 7 


+ 


20 


- 8 


+ 


52 


- 27 


+ 


50 


- 31 






+ 79 






+ 81 






.079 






.081 






285 






295 






+ .2588 






+ .4226 






- .9659 






- .9063 






+ .020 






+ .034 






- .076 






- .073 
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Table XIV (cont'd) 



II IV 





120-130 


300-310 


48 




15 


- 15 


47 


+ 


10 


- 20 


Z 


+ 


25 


- 35 


II - IV 






+ 60 


d0 






.060 


A 






305 


cos 






+ .573^ 


sin 






- .8192 


d| 






+ .034 


dn 






- .049 


46 


+ 


5 


- 25 


45 


- 


5 


- 25 


44 


- 


5 


- 20 


43 




0 


- 10 


Z 


_ 


5 


- 80 


II - IV 






+ 75 


d0 






.075 


A 






305 


cos 






+ .5736 


sin 






- .8192 


dg 






+ .043 


dTi 






- .061 


42 


+ 


7 


- 5 


41 


+ 


13 


- 8 


4o 




17 


- 14 


39 


+ 


20 


- 8 


Z 


+ 


57 


- 35 


II - IV 






+ 92 


d0 






.092 


A 






305 


cos 






+ .5736 


sin 






- .8192 


d| 






+ .053 


dTJ 






- .075 



II 


IV 


II 


130-140 


310-320 


140-150 


+ 15 


- 20 


+ 10 


+ 15 


- 15 


+ 15 


+ 30 


- 35 


+ 25 



+ 65 



.065 

315 

+ .7071 

- .7071 

+ .046 

- .046 

+ 10 - 15 +15 

+ 5 - 15 +15 

0 - 5 +10 

0 - 5 +10 

+ 15 - 4 o +45 

+ 55 

.055 
315 
+ .7071 

- .7071 

+ .039 

- .039 



+ 


5 


- 5 


+ 


8 


+ 


13 


- 10 


+ 


10 


+ 


16 


- 14 


4* 


10 


+ 


18 


- 5 


+ 


10 




52 


- 34 


4* 


38 



+ 86 



. 086 

315 
+ .7071 

- .7071 
+ . 061 

- .061 



IV 

320-330 



- 10 
- 10 

- 20 
+ 45 

.045 

325 
+ .8192 

- .5736 
+ .037 

- .026 

- 15 

- 10 
- 10 
- 10 

- 45 
+ 90 

.090 

325 

+ .8192 

- .5736 
+ .074 

- .052 

- 15 

- 20 
- 10 

- 8 

- 53 

+ 91 
.091 
325 
+ .8192 

- .5736 
+ .075 

- .052 
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Table XIV (cont'd) 







II 




IV 


II 






IV 


II 




IV 






150- 


-l6c 


330-340 


160-170 


340-350 


170- 


180 


350-360 


48 






0 




0 


. 


5 




- 10 




10 


- 20 


47 




+ 


10 




- 5 


+ 


10 




- 15 


+ 


10 


- 15 


£ 




+ 


10 




- 5 


+ 


5 




- 25 




0 


- 35 


II - 


IV 








+ 15 








+ 30 






+ 35 


d$ 










.015 








.070 






.035 


A 










335 








345 






355 


cos 








+ 


.9063 






+ 


.9659 






+ .9962 


sin 








- 


.4226 






- 


.2588 






- .0872 


di 








•f 


.014 






+ 


.029 






+ .035 


dTi 








- 


.006 






- 


.008 






- .003 


46 




+ 


15 




- 5 


+ 


15 




- 10 


+ 


25 


- 5 


45 




+ 


15 




- 10 


+ 


15 




- 10 


+ 


25 


0 


44 




+ 


10 




- 20 


+ 


10 




- 10 


+ 


15 


+ 10 


43 




+ 


5 




- 20 


+ 


5 




0 


+ 


10 


+ 15 


£ 




+ 


45 




- 55 


+ 


45 




- 30 


+ 


75 


+ 20 


II - 


IV 








+ 100 








+ 75 






+ 55 


d0 










.100 








.075 






.055 


A 










335 








345 






355 


cos 










.9063 






+ 


.9659 






+ .9962 


sin 








- 


.4226 






- 


.2588 






- .0872 


d5 








+ 


.091 






+ 


.072 






+ .055 


dT) 








- 


.042 






- 


.019 






- .005 


42 




+ 


8 




- 20 


+ 


7 




- 10 


+ 


9 


+ 12 


4 l 




+ 


10 




- 15 


+ 


10 




- 13 


+ 


12 


+ 5 


4 o 




+ 


5 




- 8 


+ 


10 




- 5 


+ 


14 


+ 4 


39 




+ 


4 




- 8 


+ 


3 




2 


+ 


10 


+ 5 


£ 




+ 


27 




- 51 


+ 


30 




- 30 


+ 


45 


+ 26 


II - 


IV 








+ 78 








+ 60 






+ 19 


de 










.078 








.060 






.019 


A 










335 








345 






355 


cos 








+ 


.9063 








.9659 






+ .9962 


sin 








- 


.4226 






- 


.2588 






- .0872 


dg 








+ 


.071 






+ 


.058 






+ .019 


dT) 








- 


.033 






- 


.016 






- .002 



r 
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Table XIV (conVd) 







I 


III 


I 




III 


I 




III 




0 - 


10 


180-190 


10 - 


20 


190-200 


20 - 


30 


200-210 


38 


+ 


17 


+ 9 


+ 


29 


+ 6 


+ 


24 


+ 3 


37 


+ 


l 6 


+ 9 


+ 


25 


+ 7 


+ 


20 


+ 5 


36 


+ 


18 


+ 1 


+ 


24 


+ 5 


+ 


19 


+ 9 


35 


+ 


21 


0 


+ 


19 


+ 1 


+ 


18 


+ 6 


£ 


+ 


72 


+ 19 


+ 


97 


+ 18 


+ 


81 


+ 23 


I - III 






+ 5 ? 






+ 79 






+ 58 


de 






.053 






.079 






.058 


A 






185 






195 






205 


cos 






- .9962 






- .9659 






- .9063 


sin 






- .0872 






- .2588 






- .4226 








- .053 






- .076 






- .053 


dTi 






- .005 






- .020 






- .024 



34 


+ 


16 


0 


+ 


19 


J 


+ 


18 


+ 2 


33 


+ 


12 


0 


+ 


16 


0 


+ 


17 


+ 2 


32 


+ 


10 


+ 2 


+ 


13 


+ 1 


+ 


16 


+ 3 


31 


+ 


8 


+ 4 


+ 


12 


+ 3 


+ 


16 


+ 2 


£ 


+ 


46 


+ 6 


+ 


60 


+ 4 


+ 


67 


+ 9 


I - III 






+ 40 






+ 56 






+ 58 


de 






.040 






.056 






.058 


A 






185 






195 






205 


cos 






- .9962 






- .9659 






- .9063 


sin 






- .0872 






- .2588 






- .4226 








- .o 4 o 






- .054 






- .053 


dTi 






- .003 






- .014 






- .025 



30 

29 

28 

27 



+ 8 +5 

+ 7 +5 

+ 6 +4 

+ 5 +4 



+ 10 +4 

+ 9 +4 

+ 7 +3 

+ 6 +3 



+ 15 +3 

+ 11 +3 

+ 9 +3 

+ 7 +3 



£ 


+ 26 +18 


+32 +14 


+ 42 +12 


I - III 


+ 8 


+ 18 


+ 30 


de 


.008 


.018 


.030 


A 


185 


195 




cos 


- .9962 


- .9659 


- .9063 


sin 


- .0872 


- .2588 


- .4226 


d| 


- .008 


- .017 


- .027 


dTi 


- .001 


- -.005 


- .013 






!• 



I 



I I 



ii 





I 



► 




t 

I 



I 
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Table XIV (cont'd) 





I 




III 


I 




III 


I 




III 




30 - 


4 o 


210-220 


4 o - 


50 


220-230 


50 - 


60 


230-240 


38 


+ 


29 


+ 7 


+ 


25 


+ 8 


+ 


16 


+ 5 


37 


+ 


22 


+ 5 


•f 


18 


+ 3 


+ 


15 


+ 1 


36 


+ 


18 


+ 5 


+ 


17 


+ 2 


+ 


15 


+ 1 


35 


+ 


17 


+ 6 


+ 


15 


+ 5 


+ 


14 


+ 5 


S 


+ 


86 


+ 23 


+ 


75 


+ 18 


+ 


60 


+ 12 


I - III 






+ 63 






+ 57 






+ 48 


d0 






.063 






.057 






.048 


A 






215 






225 






235 


cos 






- .8192 






- .7071 






- .5736 


sin 






- .5736 






- .7071 






- .8192 








- .052 






- .o 4 o 






- .028 


dTi 






- .036 






- .o 4 o 






- .039 


34 


+ 


15 


+ 6 


+ 


14 


+ 7 


+ 


13 


+ 8 


33 


+ 


15 


+ 6 


+ 


l 4 


+ 8 


+ 


14 


+ 10 


32 


+ 


17 


+ 4 


+ 


18 


+ i 


+ 


19 


+ 4 


31 


+ 


19 


+ 2 


+ 


20 


+ 1 


+ 


20 


0 


S 


+ 


66 


+ 18 


+ 


66 


+ 17 


+ 


66 


+ 22 


I - III 






+ 48 






+ 49 






+ 44 


d 0 






.048 






.049 






.o 44 


A 






215 






225 






235 


cos 






- .8192 






- .7071 






- .5736 


sin 






- .5736 






- .7071 






- .8192 








- .039 






- .035 






- .025 


dTi 






- .028 






- .035 






- .036 



30 

29 

28 

27 



+ 

+ 

+ 

+ 



l6 

12 

9 

7 

44 



+ 2 
+ 2 
+ 3 

+ 2 



+ 

+ 

+ 

+ 



l6 

12 

9 

7 

44 



+ 1 

+ 2 

+ 2 

+ 2 



+ 

+ 

+ 

+ 



17 

13 

8 

6 

44 



+ 1 

+ 2 

+ 2 

+ 2 



I - III 


+ 35 


+ 37 


+ 37 


d0 


.035 


.037 


-.037 


A 


215 


225 


235 


cos 


- .8192 


- .7071 


- .5736 


sin 


- .5736 


- .7071 


- .8192 




- .029 


- .026 


- .021 


dTi 


- .020 


- .026 


- .030 
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Table XIV (cont'd) 





I 




III 


I 




III 


I 




III 




6 o - 


70 


240-250 


TO . 


80 


250-260 


80 . 


90 


260-270 


38 


+ 


11 


0 


+ 


10 


0 


+ 


10 


+ 1 


37 


+ 


13 


0 


+ 


12 


0 


+ 


13 


+ 1 


36 


+ 


13 


+ 1 


+ 


11 


+ 2 


+ 


14 


+ 3 


35 


+ 


12 


+ 5 


+ 


10 


+ 5 


+ 


15 


+ 6 


£ 


+ 


49 


+ 6 




43 


+ 7 


+ 


52 


+ 11 


I - III 






+ 43 






+ 36 






+ 4 l 


de 






.043 






.036 






.C 4 l 


A 






245 






255 






265 


cos 






- .4226 






- .2588 






- .0872 


sin 






- .9063 






- .9659 






- .9962 








- .018 






- .009 






- .004 


dTi 






- .039 






- .035 






- . 04 l 


3^ 


+ 


12 


+ 9 


+ 


10 


+ 9 


+ 


15 


+ 6 


33 


+ 


14 


+ 10 


+ 


14 


’ 7 


+ 


15 


+ 3 


32 


+ 


19 


+ 3 


+ 


20 


+ 1 


+ 


18 


0 


31 


+ 


20 


0 


+ 


19 


0 


+ 


18 


0 


£ 


+ 


65 


+ 22 


+ 


63 


+ 17 


+ 


66 


+ 9 


I - III 






+ h3 






+ 47 






+ 57 


d0 






.043 






.0^1 






.057 


A 






245 






255 






265 


cos 






- .4226 






- .2588 






- .0872 


sin 






- .9063 






- .9659 






- .9962 


U 






- .018 






- .012 






- .005 


dn 






- .039 






- .045 






- .057 



30 


+ 


16 


+ 


1 


+ 


15 


0 


+ 


14 


0 


29 


+ 


12 


+ 


1 


+ 


10 


0 


+ 


10 


0 


28 


+ 


8 




0 


+ 


7 


0 


+ 


7 


0 


27 




6 




0 


+ 


5 


0 


+ 


5 


0 


£ 


+ 


42 


+ 


2 


+ 


37 


0 


+ 


36 


0 



I - III 


+ 40 


+ 37 


+ 36 


de 


.o 4 o 


.037 


.036 


A 


245 


255 


265 


cos 


- .4226 


- .2588 


- .0872 


sin 


- .9063 


- .9659 


- .9962 


d| 


- .017 


- .010 


- .003 


dT) 


- .036 


- .036 


- .036 
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Table XIV (conf d) 



II 

90-100 



IV II 

270-280 100-110 



IV II IV 

280-290 110-120 290-300 



38 


+ 


15 


37 


+ 


15 


36 


+ 


17 


35 


+ 


18 






65 



II - IV 

de 

A 

cos 

sin 

dn 



34 


+ 


20 


33 


+ 


19 


32 


+ 


18 


31 


+ 


16 




+ 


73 



II - IV 

de 



cos 

sin 

d| 

dr \ 



30 


+ 


13 


29 


+ 


10 


28 


+ 


7 


27 


+ 


5 




+ 


35 



II - IV 

de 



cos 

sin 

dri 



3 +20 

- 1 +20 

+ 2 +20 

+ 6 +20 

+ 4 +80 

+ 61 
.061 

275 
+ .0872 

- .9962 
+ .006 

- .061 

+ 5 +20 

+ 3 +20 

+ 1 +18 

0 +15 

+ 9 +73 

+ 64 
.064 
275 
+ .0872 

- .9962 
+ .006 

- .o 64 

0 +12 

0 +10 

0 +7 

0 +5 

0 +34 

+ 35 
.035 

275 
+ .0872 

- .9962 
+ .003 

- .035 





6 


+ 


20 


- 


5 


+ 


20 


- 


2 


+ 


20 


+ 


5 


+ 


20 




8 


+ 


80 


+ 


88 







.088 
285 
+ .2588 

- .9659 
+ .023 

- .085 



+ 


5 


+ 


20 


+ 


4 


+ 


20 


+ 


1 


+ 


18 




0 


+ 


15 


+ 


10 


+ 


73 


+ 


63 







.063 
285 
+ .2588 

- .9659 

+ .016 

- .061 

0 +12 

0 +9 

0 +6 

0 +5 

0 +32 

+ 34 
.034 
285 
+ .2588 

- .9659 

+ .009 

- .033 



- 5 

- 4 

0 

+ 5 

- 4 

+ 84 
.084 

295 
+ .4226 

- .9063 

+ .035 

- .076 

+ 5 

+ 7 

+ 2 
0 

+ 14 
+ 59 
.059 

295 

+ .4226 

- .9063 
+ .025 

- .054 

0 

0 

0 

0 

0 

+ 32 
.032 
295 
+ .4226 

- .9063 

+ .014 

- .029 
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38 

37 

36 

35 

z 

II " 

d0 

A 

cos 

sin 

d| 

dll 

34 

33 

32 

31 

Z 

II - 

d0 

A 

cos 

sin 

d| 

dT) 

30 

29 

28 

27 

z 

II - 

d0 

A 

cos 

sin 

d| 

dn 



Table XIV (cont'd) 



IV 



II 


IV 


II 


120-130 


300-310 


130-140 


+ 20 


0 


+ 19 


+ 20 


- 3 


+ 19 


+ 20 


- 1 


+ 19 


+ 20 


+ 2 


+ 20 


+ 80 


OJ CVI 
00 

1 + 


+ 77 



.082 

305 
+ .5736 

- .8192 
+ .047 

- .067 



IV 


II 


IV 


310-320 


140-150 


320-330 


- 5 


+ 10 


- 8 


- 6 


+ 15 


- 6 


2 


+ 16 


- 3 


+ 1 


+ 18 


0 


- 12 


+ 59 


- 17 


+ 89 




+ 76 


.082 




.076 


315 




325 


+ .7071 




+ .8192 


- .7071 




- .5736 


+ .063 




+ .062 


- .063 




- .o44 



+ -20 
+ 19 
+ 17 

+ 14 



+ 5 

+ 8 
+ 2 
0 



+ 20 
+ 17 
+ 14 
+ 12 



+ 3 

+ 3 

+ 2 
0 



+ 19 
+ 15 
+ 13 
+ 11 



+ 2 

+ 2 

0 
0 



+ 


70 




+ 15 


+ 


63 


+ 8 


+ 


58 




+ 4 








+ 55 






+ 55 








+ 54 








.055 






.055 








.054 








305 






315 








325 






+ 


.5736 






+ .7071 






+ 


.8192 






- 


.8192 






- .7071 






- 


.5736 






+ 


.032 






+ .039 






+ 


.044 






- 


.045 






- .039 






- 


.031 


+ 


11 




0 


+ 


11 


0 


+ 


10 




0 


+ 


9 




0 


+ 


9 


0 


+ 


8 




0 


+ 


7 




0 


+ 


7 


0 


+ 


6 




0 


+ 


5 




0 


+ 


5 


0 


+ 


5 




0 


+ 


32 




0 


+ 


32 


0 


+ 


29 




0 








+ 32 






+ 32 








+ 29 








.032 






.032 








.029 








305 






315 








325 






+ 


.5736 






+ .7071 






+ 


.8192 






- 


.8192 






- .7071 






- 


.5736 






+ 


.018 






't. .023 








.024 






- 


.026 






- .023 






- 


.017 
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Table XIV (cont'd) 



38 

37 

36 

35 

E 

II - IV 

de 

A 

cos 

sin 

«i| 

dTi 

3^^ 

33 

32 

31 

E 

II - IV 

de 

A 

cos 

sin 

d^ 

dTi 

30 

29 

28 

27 

E 

II - IV 

de 

A 

cos 

sin 

d| 

dn 



II 

150-160 



+ k 

+ 7 

+ 10 
+ 11 

+ 32 



+ 12 
+ 15 

+ 13 

+ 10 

+ 50 



+ 9 
+ 8 
+ 6 
+ 5 

+ 28 



IV 

330-340 



- 2 

- 4 

- 3 

- 1 

- 10 
+ 42 
.042 
335 
+ .9063 

- .4226 
+ .038 

- .018 

+ 1 
+ 1 
+ 1 
0 

+ 3 

+ 47 

.047 
335 
+ .9063 

- .4226 

+ .043 

- .020 

0 

0 

0 

0 

0 

+ 28 
.028 
335 
+ .9063 

- .4226 
+ .025 

- .012 



II 

160-170 



+ 1 
+ 3 

+ 4 

+ 1 

+ 9 



+ 5 

+ 10 

+ 9 
+ 8 

+ 32 



+ 8 
+ 7 
+ 5 
+ 5 

+ 25 



IV 

340-350 



0 

0 

+ 2 

+ 2 

+ 4 

+ 5 

.005 
345 
+ .9659 

- .2588 
+ .005 

- .001 

+ 3 

+ 3 

+ 2 
+ 2 

+ 10 
+ 22 
.022 
345 
+ .9659 

- .2588 

+ .021 

- .006 

+ 2 

+ 3 

+ 1 
+ 1 

+ 7 

+ 18 
.018 

345 
+ .9659 

- .2588 
+ .017 

- .005 



II 

170-180 



+ 5 

+ 4 

+ 1 
0 

+ 10 

- 18 
.018 

175 

- .9962 
+ .0872 

- .018 

+ .002 

+ 1 

+ 3 

+ 6 

+ 7 

+ 17 

- 7 

.007 

175 

- .9962 

+ .0872 

- .007 

+ .001 

+ 6 
+ 5 

+ 5 

+ 4 

+ 20 



IV 

350-360 



+ 5 
+ 5 
+ 8 
+ 10 

+ 28 



+ 8 
+ 6 
+ 5 
+ 5 

+ 24 



+ 5 

+ 5 

+ 5 

+ 3 

+ 18 
+ 2 
.002 
355 
+ .9962 

- .0872 

+ .002 

- .000 
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Table XV 


Summation 01" 


jiffectcj ol 


Inner Circle 










and Zones 


27-48 for Lai ska 




Zone 

jroup 


Group Sum 

de" 


Accumulative 

d|" 


Group Sum 
dT)" 


Accumulative 

dT)" 


Outer 
Radius 
(km. ) 


Inner 


+ 


0.010 


+ 0.010 


- 0.268 


- 0.268 


10.15 


27 - 30 


- 


0.023 


- 0.013 


- 0.383 


- 0.651 


20.09 


31 - 34 


- 


0.062 


- o’. 075 


- 0.601 


- 1.252 


39.67 


35 - 38 


- 


0.072 


- 0.147 


- 0.692 


- 1.944 


77.97 


39 - 42 


+ 


0 .l 64 


+ 0.017 


- 0.667 


- 2.611 


151.9 


43 - 46 


+ 


0.468 


+ 0.485 


- 0.763 


- 3.374 


291.2 


47-48 


+ 


0.060 


+ 0.545 


- 0.866 


- 4.240 


399.0 


2 1 

g 


= + 


0.545 


T) = 
g 


- 4 . 24 o 











# 










A/6, /5 - E/^^SCT OAJ D£^l.£C,TiaNS OP’ / NC /SEAS/a/<d 
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7»^3 Computation for Wonderfontein, Union of South Africa. 



As mentioned in section 7*^1> both Bouguer anomaly maps and 
topographic maps we 3 re required for this station. This situation 
prevailed from the outer zone, no. 48, through zone no. 31* The area 
within the circle described by that zone's radius, 20.09 km., was 
covered by a larger scale free -air anomaly map prepared by the writer. 

The inner circle for this station had a radius r^ = 5* 125 km. 
Figure l6, below is a reproduction of the inner circle. Using the 
values shown, we obtain; - 



d|" = 0.0525 (Agjj - Ag^) 

= 0.0525 ( 99.0 - 88.5) = 0.0525 (10.5) = + 0'.'551 
dn" = 0.0525 (Ag^ . Ag^) 

= 0.0525 (88.2 - 90.2) = 0.0525 (-2.0) =-0'.'105 



99.0 




Fig. 16 - Inner Circle - Wonderfontein. 
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The effects of both the inner circle and the circle-ring zones 
are accumulated in Table XVI. Curves of the effect on deflections of 
increasing gravity field radius are shown in fig. 17* 

7.^^ Comments on the Results of Deflection Computation. 

The results shown in Tables XV and XVI and in Figures 15 and 
a'', are the computed effects of the gravity field aroxmd each point 
out to a radius of 399* 0 km. ^ The effect of the rest of the world, 
beyond 399«0 km. should not be ignored in practice. Time and avail- 
ability of data do not permit careful computation of this effect for 
purposes of this thesis. However, an idea of its magnitude can be 
gotten from publications containing resvilts of investigations into 
this problem by various authors. For instance, Kaula (19)> (20), 
determined the RMS uncertainties in the total deflection as computed 
from the Vening Meinesz formulas, first assming perfect knowledge of 
gravity within given radii and none beyond, and second, ideal 
distribution of given numbers of gravity stations out to given radii. 
Taking Kaula' s figures for the first case, assuming perfect knowledge 
within a radius of 200 km. and none beyond, the RMS uncertainty is 
4Vl; for 500 km. it is 3'-'3> for 1000 km. it is 2'.'6; etc. It should 
be noted that Kaula' s computations were based on information provided 
by Hirvonen (39)- 

Interpolating for 399 4m. radius, an RMS uncertainty of + 3'*6 
is obtained if Kaula' s information is assumed to apply. The probable 
error is then + 0 . 67^5 (3'.6) = + 2','U. 
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Table XVI - Summation of Effects of Inner Circle 
and Zones 23 - ^ for Wonderfontein. 



Zone 


Group Sum 


Accumulative 


Group 




d|" 


d|" 


Inner 


+ 


0.551 


+ 0.551 


23 - 


26 


+ 


0.515 


+ 1.066 


27 - 


30 


+ 


0.905 


+ 1.971 


31 - 




+ 


0.128 


+ 2.099 


35 - 


38 


+ 


0.175 


+ 2.274 


39 - 


h 2 


- 


1.884 


+ 0.390 


i+3 - 


k 6 


- 


1.919 


- 1.529 


U? - 


48 


- 


0.75^^ 


- 2.283 


Z 




= - 


2.283 





Group Sum 
dTi" 


Accumulative 

dT)" 


Outer 
Radius 
(km. ) 


- 0.105 


- 0.105 


5.125 


+ 0.013 


- 0.092 


10.15 


+ 0.212 


+ 0.120 


20.09 


+ 1.099 


+ 1.219 


39.67 


+ 1.448 


+ 2.667 


77.97 


+ 1.139 


+ 3.806 


151.9 


+ 1.222 


+ 5.028 


291.2 


- 0.493 
= + 4.535 


+ 4.535 


399.0 



f 




'-1 
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Estimates of the probable errors of the gravimetric deflections 
as computed from the data given within the field of radius 399*0 km. 
must be made, since perfect knowledge of gravity within the field can 
not be assumed. The method to be used for this estimate is that of 
the U.S.C. & G.S. and is described in Rice's paper on the l6 points 
(26). According to this method, the probable error in either component 
i or Ti, is determined from the following expression; - 

= 18 z (o'.'oo;l 

where E^ = Probable error in determining the mean anomaly 

in any compartment in a sector in mgals. 

E^ = Probable error in | or t) in seconds 

In making use of the above equation for this problem, the 
probable error was estimated for each zone group computed. The 
accumulated results, combined with the effect of the innef~ circle, each 
component of which is assumed to be no more certain than its value, and 
with the probable error for the area beyond the radius 399* 0 km. , give 
the following (see tables XVII and XVIIl):- 

For Total Deflection 9 

Laiska \ “ i 2 V 42 

Wonderfontein = + 2 '.' 47 . 

Thus it appears that for computations based upon the gravity 
data given, and within a field of radius 399*0 km., the probable eirors 
for the gravity field are almost negligible compared with those re- 
sulting from the neglect of gravity beyond that radius. 
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Table XVII - Computation of Probable Errors in Gravimetric 

Deflections as computed for the Gravity- 
Field within Radius 399«Q tan. 



(1) 


(2) 


(3) 




(5) 


(6) 


(7) 


• 






E 

Agg 

(mgal) 




E 

r 


0.001 E 
r 


(O.OOl E f 
r 


2(O.U01E_)‘ 
' r 


Zones 


No. of 
Zones 


(mgal) 


(rigal) 






for Group 










i[(3)t(‘>)'l 0.001(5) 


(6)^ 


(2) (7) 


LAISKA 
















27-30 


4 






’ 1 


0.001 


0.000001 


0. 000004 


31-3^ 


4 






1 


.001 


.000001 


. 000004 


35-38 


4 






1 


.001 


. 000001 


. 000004 


39-^2 


4 






5 


.005 


. 000025 


.000100 


43-46 


4 






5 


.005 


.000025 


. 000100 


47-48 


2 






5 


.005 


. 000025 


.000050 



WONDERFONTEIN 
23-26 4 1 


r 3 


2 


0.002 


0.000004 


0. 000016 


27-30 


4 


1 


i 3 


2 


.002 


. 000004 


.000016 


31-34 


4 


1 


L 3 


2 


.002 


. 000004 


.000016 


35-38 


4 


5 


2{ 7 


6 


.006 


. 000036 


. 000144 


39-^2 


4 


5 


fl7 


11 


.011 


.000121 


. 000484 


43-46 


4 


5 


3 j 17 


11 


.011 


. 000121 


.000484 


47-48 


2 


5 


LI7 


11 


.011 


. 000121 


. 000242 


1 100 
2 200 


ft. 

ft. 


@ 0.03378 mgal / ft. 
@ 0.03378 mgal / ft. 









3 500 ft. @ 0.03378 mgal / ft. 



Inner 


Circles ;- 










Laiska 




Wonderfontein 




; E. = d| + OVOlO 

2 

®t 


0.000100 


+ 0.551 


0. 303601 




: E^ = dTi - 0.268 

^2 
h 


0.071824 


- 0V105 


0.011025 




m 2 .,2 , ,^2 

: E^ = d| + dT) 


0.071924 




0.31^4626 



Table XVII (cont*d) 



l8 2 ( 0.001 Er) j ~ p.e. in either component 

(1 or Ti) of the deflection 



(9) 

Zones 



( 10 ) 
18 (8) 



( 11 ) 



Accumulated 
„ 2 



( 12 ) 

Accumulated 

3 — 

^ (11) 



LAISKA 








2'f-30 


0.000072 


0.000072 


+ o'.'oo8 


31-34 


.000072 


.000144 


.012 


35-38 


.000072 


.000216 


.015 


39-^'-2 


.001800 


.002016 


.045 


43-46 


.001800 


.003816 


.062 


47-48 


.000900 


.004716 


.068 



For Total Deflection 0 : 



E, 



2(0.004716) 
0. 009432 
o"o97 



WOKDERFONTEIN 



23-26 

27-30 

31-34 

35-38 

39-42 

43-46 

47-48 



0. 000088 
. 000288 
. 000288 
.002592 
.008712 
.008712 
. 004356 



For Total Deflection 0 : E. 



E^ 



0.000288 + 0'.'017 

.000576 .024 

.000864 .029 

.003456 .059 

.012168 .110 

.020880 .145 

.025236 .159 

= 2(0.025236) 

= 0.050472 
= 0'.'225 
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Table XVIII - Probable 


Errors 


in Gravimetric 


Deflections 1 




for 


Combined Zones 




Zones 




Accumulated 

0 


Accumulated 








\ 


LAISKA 








Inner Circle = 10.25 km. ) 


0.0719 


0.0719 


j'.'268 


27-48 ( 10.25 - 399.0 km.) 


0. 0094 


0.0813 


0"285 


Beyond 399*0 km. ; ( 2 '.' 4 )^ 


5.76 


5.8413 


2'.' 42 


WONDERFOOTEIN 








Inner Circle (r^ = 5*125 km.) 


0 . 31^*6 


0.3146 


0'.'56i 


23-26 (5.125 - 399.0 km.) 


0.0505 


0. 3651 


0'.'6o4 


Beyond 399-0 km. : ( 2 '.' 4 )^ 


5.76 


6.1251 


2 '.' 47 



I 



I 



i 
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7«5 Approximate Correction to Finland - South Africa Arc from 

Gravimetric Deflections 

Making use of the values obtained for the gravimetric deflect- 
ions of Laiska and Wonderfontein in section 7*^# disregarding the un- 
certainty arising from the fact that computation was not carried out 
beyond 399 tan., a correction to the length of the meridian arc between 
the points can be made. From this, a rovjgh idea of the corresponding 
correction to the equatorial radius can be learned. 

Geodetic and astronomic data for the two points in question 
was obtained from the Army Map Service, and was part of the data used 
in that agency's computation of "A New Determination of the Figure of 
the Earth from Arcs" (See Chovitz and Fischer (4)). Unfortunately the 
data carries a security classification for values having accuracy 
better than one minute of arc. Because of this, the actual computation 
figures cannot be published here. However, the results given reflect 
the data furnished to the writer by AMS without violating security 
regulations. 

Using latitude fxanctions tables (either (34) or (36)) for the 
International Ellipsoid, upon which the given data was based, the 
length of the meridian arc between Laiska and Wonderiontein for the 
geodetic latitudes and for the latitudes corrected by the free-air 
gravimetric and astro-geodetic deflections, were computed (See figure 
l8). In effect, the two figures thus obtained were (l) the measured 
length of the arc, based on the International Ellipsoid, and (2) the 
corrected length of the arc based on astro-gravimetric data. The 
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following figures give pertinent details 



Station 



Latitude 



Meridian Arc 



Laiska 

Wonderfontein 






M 



cp = 64° 03' ( )" N 

cp = 25° 48' ( )" S 

Ziscp = 89° 51’ ( )" 

(M^ + M^) = 9,961,4^8 m. 



Laiska 

Wonderf onte in 



cp = 64° 03* ( )" N 

s ^ 

cp^ = 25° 48' ( )" S 
Acp = 89° 51' ( )" 

(\ ♦ \)g 



(«L)g 



9,561,5311 m. 



where Mj^ and refer to the meridian distance from the equator to the 

point in question, cp is the geodetic latitude, and cp the graviraetric- 

S 

ally corrected latitude. Thus cp = cp+(9'"'?)“6* 

S S 

From the above figures, we see that the total correction for 
the arc is + 86 meters. Assuming the arc to be a meridian arc of a 
sphere the size of the earth, we can thus obtain a rough correction to 
the sphere's radius 



5R = 



£_ 

zkp 



86 



89?85 



= + 55 meters 6a 



Admittedly, this is only a rough approximation, but is close 
enough to the exact value for the data used, for illustrative purposes. 
Besides this, the accuracy of the gravimetric deflections is of doubt- 



€•• - 
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ful value, because of the limited gravity field used. 

Any person who is lucky enough to be in possession of data 
such as is used here, but for many points along an arc of this length, 
would proceed along the lines of section 7*2, and would doubtless 
obtain results of some significance. 



NO/tZTH 

/=>Oi^E 




Pig. 18 



The Finland - South Africa Arc 
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Chapter 8 - Closing Remarks 

It would be presumptuous on the part of the writer to title 
this chapter "Conclusions", because there can hardly be any conclusion 
to something that has not, in a manner of speaking, even begun. The 
lack of readily available data for pursuing a project of such 
magnitude as determining the best dimensions for the earth ellipsoid, 
prevents a significant contribution to the science, even if the ability 
to put such data to proper use were present. 

What conclusions can be made, must be based on what was al- 
ready known beforehand, and indeed could not be made without such 
knowledge, because of the nature and amount of the data used. The 
first two problems solved in Chapter 7 are what are referred to here. 

The use of only astro-geodetic deflections of the vertical is admitted- 
ly an incomplete beginning to a solution, so could hardly be expected 
to yield as good results as with the use of the gravimetric deflections. 
It is unfortunate that gravimetric deflections were not available for 
all the points used for the arc of the 98th meridian. A direct com- 
parison of methods could then have been made. It would have been 
almost as meaningless to have computed the size of the earth from only 
the astro-geodetic deflections of the 16 points in the South Central 
U. S., because of the small area covered. As it was, the computation 
by the astro-gravimetric method in the case of the 16 points, probably 
a/ed much of its "success" to the fact that the data given, as 
prepared by Mr. Rice, vms carefully computed by him and fairly 
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accurate. 

It is hoped that the material presented in this thesis might 
provide, by its systematic presentation of the procedure to be 
followed in problems of this nature, the means whereby any interested 
reader might pursue the problem further, especially when in the future 
more gravimetric deflections of the vertical can be expected to be 
available. The advent of the high speed electronic computer in 
geodetic use will hasten tha-f day. The manual method of computing 
gravimetric deflections is so laborious that one person, in the vrork 
of preparation of a thesis, could not be expected to produce the 
required computations in adequate numbers and for those stations in 
locations favorable to solution of the problem of determining the size 
of the earth ellipsoid. It is believed that the example deflection 
computations included in this thesis demonstrate this clearly. 

It is regretted that more them mere mention of two important 
methods used for solution of the problem at hand could not have been 
made in this thesis. One is the determination of the figure of the 
earth from geoid heights, and the other is the determination of the 
flattening from observations of the artificial earth satellites. 

Time, more than space, prevented their inclusion, especially since the 
writer tends to be somewhat verbose in his explanations. 
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